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ABSTRACT 


This  thesis  explores  the  use  of  a  pointed  seeker  lens 
designed  using  a  spherically  symmetric  gradient  refractive 
index  (GRIN) .  The  design  helps  to  solve  the  current  design 
conflict  between  optical  quality  and  aerodynamic  drag  in¬ 
herent  in  hemispherical  seeker  lenses.  Equations  for  lens 
design  and  the  evaluation  of  off-axis  lens  performance  have 
been  developed  for  both  a  homogeneous  version  and  a  GRIN 
version  of  the  pointed  seeker  lens.  The  homogeneous  lens 
is  used  as  a  comparison  and  a  check  for  the  GRIN  lens.  A 
FORTRAN  program  (GISL)  has  been  written  and  employed  to 
evaluate  and  compare  both  the  homogeneous  lens  and  many  dif¬ 
ferent  configurations  of  possible  GRIN  lens  designs.  Results 
indicate  that  the  GRIN  lens  has  highly  superior  off-axis  imag¬ 
ing  performance  as  compared  to  the  homogeneous  lens.  The  best 
results  were  obtained  for  the  GRIN  lens  with  a  fifty  percent, 
positive,  spherically  symmetric  gradient  index  with  center  of 
symmetry  interior  to  the  lens.  Only  very  slightly  inferior 
performance  was  observed  with  a  five  percent  version  of  the 
same  lens;  such  a  lens  possibly  can  be  manufactured  today. 

GRIN  lens  performance  also  indicates  that  for  objects  off- 
axis  by  more  than  17.2  degrees  a  large  scale,  multiple  ele¬ 
ment  sensor  array  may  be  required;  with  such  a  sensor  array, 
objects  off-axis  by  more  than  37.2  degrees  may  require  mirror 
elements  to  compensate  for  image  movement. 
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Spot  Diagram  for  Grid  of  Figure  F-30  - 

Encircled  Energy  of  Figure  F-31  - 

GRIN  Lens  Shape  at  -5%,  OB  3  0.20,  a  3  9.00  — 

Grid  Plane  at  ap  3  0.3  for  Lens  of  Figure  F-33 

Spot  Diagram  for  Grid  of  Figure  F-34  - - 

Encircled  Energy  of  Figure  F-35  — - - - - - 

GRIN  Lens  Shape  at  -10%,  OB  3  0.20,  a  3  9.00  - 

Grid  Plane  at  ap  3  0.3  for  Lens  of  Figure  F-37 

Spot  Diagram  for  Grid  of  Figure  F-38  — - - 

Encircled  Energy  of  Figure  F-39  - - - 

GRIN  Lens  Shape  at  -25%,  OB  3  0.20,  a  3  9.00  - 

Grid  Plane  at  ap  3  0.3  for  Lens  of  Figure  F-41 

Spot  Diagram  for  Grid  of  Figure  F-42  — - — - 

Encircled  Energy  of  Figure  F-43  - - - - - 

GRIN  Lens  Shape  for  -50%,  OB  3  0.20,  a  3  9.00 


F-46.  Grid  Plane  at  a  *  0.3  for  Lens  of  Figure  F-45 

P 

F-47.  Spot  Diagram  for  Grid  of  Figure  F-46  - 

F-48.  Encircled  Energy  of  Figure  F-47  - — 

F-49.  GRIN  Lens  Shape  for  +5%,  OB  *  0.20,  a  =  9.00  - 

F-50.  Grid  Plane  at  a  *0.3  for  Lens  of  Figure  F-49 

P 

F-51.  Spot  Diagram  for  Grid  of  Figure  F-50  - 

F-52.  Encircled  Energy  of  Figure  F-51  - — 

F-53.  GRIN  Lens  Shape  at  4-10%,  OB  *  0.20,  a  *  9.00  - 

F-54 .  Grid  Plane  at  ap  *  0.3  for  Lens  of  Figure  F-53 

F-55 .  Spot  Diagram  for  Grid  of  Figure  F-54  - - 

F-56 .  Encircled  Energy  of  Figure  F-55  - 

F-57.  GRIN  Lens  Shape  at  4-25%,  OB  -  0.20,  a  *  9.00  - 

F-58 .  Grid  Plane  at  ap  ■  0.3  for  Lens  of  Figure  F-57 

F-59 .  Spot  Diagram  for  Grid  of  Figure  F-58  - 

F-60 .  Encircled  Energy  of  Figure  F-59  - - 

F-61.  GRIN  Lens  Shape  at  4-50%,  OB  -  0 .20 ,  a  -  9 .00  - 

F-62.  Grid  Plane  at  <*  *  0.3  for  Lens  of  Figure  F-61 

P 

F-63.  Spot  Diagram  for  Grid  of  Figure  F-62  - 

F-64.  Encircled  Energy  of  Figure  F-63  - - - - 

F-65.  GRIN  Lens  Shape  at  -5%,  OB  ■  1.00,  a  »  9.00  — 
F-66 .  Grid  Plane  at  a  *  0.3  for  Lens  of  Figure  F-65 

F-67 .  Spot  Diagram  for  Grid  of  Figure  F-66  - — 

F-68 .  Encircled  Energy  of  Figure  F-67  - - - 

F-69 .  GRIN  Lens  Shape  at  -10%,  OB  *  1.00,  a  »  9.00  - 
F-70 .  Grid  Plane  at  a  ■  0.3  for  Lens  of  Figure  F-69 
F-71.  Spot  Diagram  for  Grid  of  Figure  F-70  - — 


Encircled  Energy  of  Figure  F-71  - - - - 

GRIM  Lens  Shape  at  -25%,  OB  ■  1.00,  a  *  9.00  - 

Grid  Plane  at  ctp  *  0.3  for  Lens  of  Figure  F-73 

Spot  Diagram  for  Grid  of  Figure  F-74  - 

Encircled  Energy  of  Figure  F-75  - - — — 

GRIN  Lens  Shape  at  -50%,  OB  *  1.00,  a  *  9.00  - 

Grid  Plane  at  a  *0.3  for  Lens  of  Figure  F-77 
P 

Spot  Diagram  for  Grid  of  Figure  F-78  - 

Encircled  Energy  of  Figure  F-79  - - 

GRIN  Lens  Shape  at  +5%,  OB  *  1.00,  a  *  9.00  -— 

Grid  Plane  at  a  *  0.3  for  Lens  of  Figure  F-81 
P 

Spot  Diagram  for  Grid  of  Figure  F-82  - - 

Encircled  Energy  of  Figure  F-83  - 

GRIN  Lens  Shape  at  +10%,  OB  *  1.00,  a  *  9.00  - 
Grid  Plane  at  a  *  0.3  for  Lens  of  Figure  F-85 

tr 

Spot  Diagram  for  Grid  of  Figure  F-86  - 

Encircled  Energy  of  Figure  F-87  - 

GRIN  Lens  Shape  at  +25%,  OB  *  1.00,  a  *  9.00  — 

Grid  Plane  at  ap  *  0.3  for  Lens  of  Figure  F-89 

Spot  Diagram  for  Grid  of  Figure  F-90  - 

Encircled  Energy  of  Figure  F-91  — - - 

GRIN  Lens  Shape  at  +50%,  OB  *  1.00,  a  *  9.00  - 

Grid  Plane  at  an  ■  0.3  for  Lens  of  Figure  F-93 
P 

Spot  Diagram  for  Grid  of  Figure  F-94  - - 

Encircled  Energy  of  Figure  F-95  - — - - 

GRIN  Lens  Shape  at  -5%,  OB  *  2.00,  a  *  9.00  — 


0.3  for  Lens  of  Figure  F-97 


F-98.  Grid  Plane  at  a  * 

P 

F-99.  Spot  Diagram  for  Grid  of  Figure  f-98  - 

F-100.  Encircled  Energy  of  Figure  F-99  - — 

F-101.  GRIN  Lens  Shape  at  -10%.  OB  =  2.00,  a  =  9.00  — 

F-102 .  Grid  Plane  at  *  0.3  for  Lens  of  Figure  F-101 

F-103.  Spot  Diagram  for  Grid  of  Figure  F-102  - 

F-104.  Encircled  Energy  of  Figure  F-103  - - - 

F-105.  GRIN  Lens  Shape  at  -25%,  OB  »  2.00,  a  *  9.00  — 

F-106 .  Grid  Plane  at  a  =»0.3  for  Lens  of  Figure  F-105 

P 

F-107.  Spot  Diagram  for  Grid  of  Figure  F-106  - 

F-108.  Encircled  Energy  of  Figure  F-107  — — - - 

F-109.  GRIN  Lens  Shape  at  -50%,  OB  =  2.00,  a  *  9.00  — 

F-110 .  Grid  Plane  at  a  =0.3  for  Grid  of  Figure  F-109 

F-lll .  Spot  Diagram  for  Grid  of  Figure  F-110  — — - — 

F-112 .  Encircled  Energy  of  Figure  F-lll  - — — 

F-113.  GRIN  Lens  Shape  at  +5%,  OB  =  4.00,  a  «  9.00  - — 

F-114.  Grid  Plane  at  a  =0.3  for  Lens  of  Figure  F-113 

P 

F-115.  Spot  Diagram  for  Grid  of  Figure  F-114  - 

F-116.  Encircled  Energy  of  Figure  F-115  - 

F-117 .  GRIN  Lens  Shape  at  +10%,  OB  =  4.00,  a  =  9.00  — 

F-118.  Grid  Plane  at  a  =  0.3  for  Lens  of  Figure  F-117 

P 

F-119.  Spot  Diagram  for  Grid  of  Figure  F-118  - 

F-120 .  Encircled  Energy  of  Figure  F-119  - - - 

F-121.  GRIN  Lens  Shape  at  +25%,  OB  =  4.00,  a  =  9.00  — 
F-122.  Grid  Plane  at  =  0.3  for  Lens  of  Figure  F-121 
F-123 .  Spot  Diagram  for  Grid  of  Figure  F-122  - — 


Encircled  Energy  of  Figure  F-123  — - — — — - — 

GRIM  Lens  Shape  at  +50%,  OB  *  4.00,  a  *  9.00  — 

Grid  Plane  at  a  =  0.3  for  Lens  of  Figure  F-125 
P 

Spot  Diagram  for  Grid  of  Figure  F-126  - 

Encircled  Energy  of  Figure  F-127  — — — - - 

"Best"  GRIN  Lens  Shape  with  50%  Gradient, 

OB  »  0.05  and  a  *  9.00  in  the  F/2  Configuration 

Grid  Plane  at  ctp  »  0.0  for  Lens  of  Figure  G-l  - 

Spot  Diagram  for  Grid  of  Figure  G-2 - - 

Grid  Plane  at  ap  *  0.1  for  Lens  of  Figure  G-l  - 

Spot  Diagram  for  Grid  of  Figure  G-4 - - — 

Encircled  Energy  of  Figure  G-5  - 

Grid  Plane  at  ap  =  0.2  for  Lens  of  Figure  G-l  - 

Spot  Diagram  for  Grid  of  Figure  G-7  - 

Encircled  Energy  of  Figure  G-8  — - - — - 

Grid  Plane  at  ap  ■  0.4  for  Lens  of  Figure  G-l  - 

Spot  Diagram  for  Grid  of  Figure  G-10  - - - 

Encircled  Energy  of  Figure  G-ll  - - — - - 

Grid  Plane  at  a  *0.5  for  Lens  of  Figure  G-l  - 
P 

Spot  Diagram  for  Grid  of  Figure  G-l 3  - - - 

Encircled  Energy  of  Figure  G-l 4  - — - - — - — — 

Grid  Plane  at  a  *  0.6  for  Lens  of  Figure  G-l  - 

Sr 

Spot  Diagram  for  Grid  of  Figure  G-l 6 - - 

Encircled  Energy  of  Figure  G-17  - - - 

Grid  Plane  at  a  *  0.7  for  Lens  of  Figure  G-l  - 

Spot  Diagram  for  Grid  of  Figure  G-19  - - — - 

Encircled  Energy  of  Figure  G-20  - - — - - 


Grid  Plane  at  ap  *  0.8  for  Lens  of  Figure  G-l  — 

Spot  Diagram  for  Grid  of  Figure  G- 22  - 

Encircled  Energy  of  Figure  G- 23  - - - - 

"Best"  GRIN  Lens  Shape  with  50%  Gradient, 

OB  »  0.05,  and  a  *  9.00  in  the  F/l  Configuration 

Grid  Plane  at  a  *  0.0  for  Lens  of  Figure  H-l  — 
P 

Spot  Diagram  for  Grid  of  Figure  H-2 - - 

Grid  plane  at  a  ■  0.1  for  Lens  of  Figure  H-l  — 
P 

Spot  Diagram  for  Grid  of  Figure  H-4 - — 

Encircled  Energy  of  Figure  H-5  - - - - - - 

Grid  Plane  at  a  *  0.2  for  Lens  of  Figure  H-l  — 
P 

Spot  Diagram  for  Grid  of  Figure  H-7  - - 

Encircled  Energy  of  Figure  H-8  - - - - 

Grid  Plane  at  a  *  0.3  for  Lens  of  Figure  H-l  — 

tr 

Spot  Diagram  for  Grid  of  Figure  H-10  - 

Encircled  Energy  of  Figure  H-ll  - - - — 

Grid  plane  at  a  *  0.4  for  Lens  of  Figure  H-l  — 
P 

Spot  Diagram  for  Grid  of  Figure  H-13  - - - - 

Encircled  Energy  of  Figure  H-14  - 

Grid  Plane  at  oip  *  0.5  for  Lens  of  Figure  H-l  — 

Spot  Diagram  for  Grid  of  Figure  H-16 - - 

Encircled  Energy  of  Figure  H-17  - - - — — 

Grid  Plane  at  a  -  0.6  for  Lens  of  Figure  H-l  — 
P 

Spot  Diagram  for  Grid  of  Figure  H-19 - — — - 

Encircled  Energy  of  Figure  H-20  — - - - - 

Grid  Plane  at  ap  *  0.7  for  Lens  of  Figure  H-l  — 
Spot  Diagram  for  Grid  of  Figure  H-22 - - 
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Defined  by  Equation  44 
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po 
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by  Figure  5 
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Vector  in  the  direction 
of  refracted  ray 
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Nondimensional 

S 

51 
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Figure  5 

Nondimensional 

ST 
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see  Figure  5 

Nondimensional 

T 
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lens 
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with  respect  to  the 
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UP 
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inside  the  lens  with 
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Radians 
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Integration  variable 
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22 

General  x-coordinate 
along  lens  axis 

Nondimensional 

30 
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(tilted) 
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2 

Absolute  value  of  any 
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XD 
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x-coordinate  of  outside 
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intercept 

Nondimensional 

XI 
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x-coordinate  of  inside 

Nondimensional 
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surface  intercept 
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x-ooordinate  in  the 
image  plane  (Spot 
Diagram) 

Nondimsnsional 
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First  x-coardinate  of 
imaginary  homogeneous 
intercept  in  GRIN 
design 
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XP 

XP 
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x-ooardinate  of  the 
first  intermediate 
point  an  the  GRIN 
ray  during  iteration 
in  the  design  algorithm 
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imaginary  homogeneous 
intercept  in  GRIN 
design 

Nondimens  ional 
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XP 
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Successive  values  of 
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x-ooordinate  of  ray 
intercept  in  meridian 
plane 
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plane 
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r 
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y-coordinate  in  the 
tilted  grid  plane 
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*0 

Y0 
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y-coordinate  of  outside 
surface  skew  ray 
intercept 
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y-coordinate  of  inside  Nondimensional 
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intercept 
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y-ooordinate  of  the 
star  ray  in  the  image 
plane  (Spot  Diagram) 
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intercept  in  GRIN  de¬ 
sign  (iteration) 
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YP 
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x-coordinate  of  the 
first  intermediate 
point  on  the  GRIN  ray 
during  iteration  in  the 
design  algorithm 

Nondimensicnal 
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Successive  values  of 

y  during  iteration 
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y-ooordinate  of  the 
image  oentroid  in  the 
Spot  Diagram 

Nondimensicnal 
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Individual  y-coordinates 
of  image  plane  ray 
intercept 

Nondimensicnal 
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y-ooordinate  of  the 

Jth  ray  on  the  outside 
surface  in  the  meridian 
plane 

Nondimensicnal 
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5 

y-ooordinate  of  the  ray 
on  the  inside  surface  in 
the  meridian  plane 

Nandi  mens ional 
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General  z-coordinate 
(horizontal) 

Nondimensicnal 
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z-coordinate  in  the 
grid  plane  (tilted) 

Nondimensicnal 

ZO 
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z -coordinate  of  the 
skew  ray  on  the  out¬ 
side  surfaoe 

Nondimensicnal 

ZI 
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z -coordinate  of  the 
skew  ray  on  the  inside 
surface  (cone) 

Nondimensicnal 

■*.  *\  - 

*'”’*■* 
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zim 

ZIM 
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z -coordinate  of  the 
skew  ray  in  the  image 
plane  (Spot  Diagram) 

Nondimens ional 
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with  respect  to  the 
lens  axis 
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Total  angle  between  lens 
axis  and  the  point  of 
maximun  radius  of  the 
cone;  measured  at  the 
focal  point 

Radians 

\ 

Y 

GNMA 

29 

Nose  half-angle 
(opaque  region) 

Radians 
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Defined  by  Equation  124 

i 

AU" 

primp 

3 

Angle  between  successive 
rays  at  the  focal  point 

Radians 

i  03 

Ax 

52 

Incremental  change  in 
x;  see  Figure  5 

Ay 
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Incremental  change  in 
y;  see  Figure  5 
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EPSILCN 
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Sign  function  (±1) 
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ZETA 

134 

Angle  between  instantan¬ 
eous  GRIN  ray  direction 
and  the  lens  axis  in  the 
meridian  plane 

Radians 

n 

120 

Defined  by  Equation  123 

1 

M* 

0 

108 

Generalized  angular  co¬ 
ordinate  of  GRIN  ray  in 
the  plane  of  the  ray 

Radians 

1 

91 

1 

Generalized  angle  of 
incidence  with  respect 
to  the  local  normal 

Radians 

W* 
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' j  s' 

92 

1 

Generalized  angle  of 
refraction  with  respect 
to  the  local  normal 

Radians 
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intercept)  in  GRIN 

Radians 
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Iteration  trial  value 
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Total  angular  (SIN 
coordinate 

Radians 

THEEAP 
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Trial  values  of  8 
during  iteration 

Radians 

XN EW 
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Revised  trial  value 
of  8  during  iteration 

Radians 

54 

Angle  between  meridian 
plane  and  the  point  of 
intersection  on  the 
outside  surface 

Radians 

PI 
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Proportionality  factor 
between  the  circun- 
ference  and  the 
diameter  of  a  circle; 
3.14159... 

Radians 

RMSRAD 
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Spot  size;  RM3  radius 
of  image 

Nondimens  ionaLL 

SIQflff 
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Standard  deviation  of 
y-ooordinates  of  rays 
in  the  Spot  Diagram 
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SIGMAZ 
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z-standard  deviation 
of  rays  in  the  Spot 
Diagram 
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57 

Angle  between  skew  ray 
and  the  outside  surface 
with  respect  to  the 
surface  normal 

Radians 
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61 

Refraction  angle  be¬ 
tween  skew  ray  and  the 
outside  surface  normal 

Radians 
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90 

Angle  of  skew  ray  inter¬ 
cept  with  respect  to  the 
inside  surfaae  normal 

Radians 
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Refraction  angle  of  skew  Radians 
ray  with  respect  to  the 
inside  surface  normal 
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intercept  at  outside 
surface 


Radians 


Radians 


34 


ACKNOWLEDGEMENTS 


The  author  would  like  to  express  his  sincere  apprecia¬ 
tion  to  Distinguished  Professor  Allen  E.  Fuhs  for  his 
invaluable  guidance,  friendship,  and  time  away  from  sabbati 
cal.  Without  him  this  thesis  would  not  have  been  possible. 

The  author  would  also  like  to  thank  Dr.  Oded  Amichai 
for  his  friendship  and  his  assistance  in  many  areas  most 
notably  in  the  intricacies  of  the  IBM  3033. 


I .  INTRODUCTION 


Historically,  the  design  of  tactical  missiles  employing 
passive  or  semiactive  infrared  (IR)  seekers  has  involved  a 
difficult  compromise  between  aerodynamic  requirements  and  opti¬ 
cal  or  seeker  requirements.  Whereas  aerodynamically  the  mis¬ 
sile  nose  region  should  be  sharp  in  order  to  reduce  drag, 
optically  it  should  be  hemispherically  shaped  for  image 
quality  and  as  large  as  possible  to  increase  aperture  and 
therefore  acquisition  or  tracking  range.  Some  IR  homing  mis¬ 
siles  designed  for  very  short  range  anti-armor  missions  have 
totally  ignored  nose  drag  in  order  to  optimize  seeker  per¬ 
formance  while  other  designs  for  longer  range  missiles  re¬ 
quiring  high  cruise  velocities  and  greater  aerodynamic 
efficiency  have  used  the  blunted  ogive  as  a  compromise.  There 
have  not  been  any  IR  designs  which  have  ignored  optics  in  favor 
of  aerodynamics;  nor  has  there  been  employed  a  pointed  seeker 
lens  with  the  desired  optical  qualities. 

In  order  to  increase  the  performance  of  optically  guided 
missiles  beyond  the  current  state  of  development,  the  conflict 
between  aerodynamic  requirements  and  optical  restrictions  must 
be  resolved.  Significant  improvement  in  missile  performance 
by  increasing  thrust  is  not  likely  due  to  the  highly  advanced 
state  of  propulsion  today.  One  way  to  resolve  the  aerodynamics 
optics  problem  is  to  design  a  pointed  lens  which  has,  if  not 


imaging  quality,  enough  optical  performance  to  allow  the 
reduction  of  tracking  data. 

The  lens  must  be  in  the  general  shape  of  a  cone  or  ogive 
with  half  angle  compatible  with  the  design  free  stream  Mach 
number  in  order  to  maintain  an  attached  shock  wave.  Such  a 
spike-shaped  lens  might  also  be  used  in  the  diffuser  portion 
of  an  optically  guided  ramjet  with  nose  inlet  to  conserve 
stagnation  pressure.  A  diffuser- lens  may  benefit  by  a  semi- 
isentropic  spike  shaped  lens. 

Poor  optical  performance  due  to  the  pointed  shape  of 
sharp  lenses  has  precluded  their  use.  Gradient  Refractive 
Index  (GRIN)  materials,  however,  permit  the  lens  designer 
the  freedom  to  spatially  vary  the  lens  index  of  refraction 
to  compensate  for  a  traditionally  poor  optical  shape.  Al¬ 
though  the  use  of  GRIN  has  seen  widespread  use  in  fiber  optics 
technology,  it  has  been  used  infrequently  in  lens  applications 
until  recently.  Lens  designers  are  discovering  that  multiple- 
element  photographic  objectives  may  be  redesigned  using  a 
two-element  gradient  lens  [1].  It  should  be  noted,  however, 
that  such  GRIN  lens  have  not  yet  been  successfully  fabricated 
even  though  large  index  changes  in  glass  have  been  accomplished 
by  the  diffusion  of  doped  electropolarizable  ions  [2] .  Ex¬ 
tensive  research  is  being  conducted  in  the  creation  of  ever 
larger  and  more  precise  gradients.  A  total  change  in  refrac¬ 
tive  index  of  approximately  five  percent  is  currently 


attainable. 


At  the  Naval  Postgraduate  School,  Frazier  [3],  Terrell 
[4],  and  Amichai  [5]  have  studied  GRIN  as  applied  to  the 
sharp  lens  problem.  The  brief  introduction  by  Frazier  was 
followed  by  Terrell  who  designed  a  sharp  lens  having  a  coni¬ 
cal  outside  surface  and  a  variable  inner  surface  using  a 
homogeneous  index  (HIN)  and  then  also  briefly  touched  on  the 
GRIN  application.  Amichai  extended  Terrell's  lens  to  the 
GRIN  case  and  included  preliminary  results  from  this  thesis 
in  a  computer  routine  intended  as  a  framework  for  lens 
optimization  by  following  researchers. 

This  thesis  investigates  the  design  and  off-axis  per¬ 
formance  (skew  rays)  of  a  GRIN  seeker  lens  having  a  variable 
outer  surface  and  a  fixed  right  circular  cone  as  an  inside 
surface.  The  variable  outside  surface  is  determined  by  the 
character  of  the  spherical  gradient  employed  and  varies  from 
a  pseudo  ogive  to  an  approximate  isentropic  spike  as  shown  in 
Figure  1.  First,  the  theory  of  the  HIN  lens  is  developed, 
followed  by  the  design  of  a  homogeneous  lens  intended  both 
as  a  comparison  and  a  check  for  the  GRIN  lens  theory  and  the 
design  which  follows.  Lens  performance  parameters  are  dis¬ 
cussed  and  results  presented  for  both  the  HIN  and  GRIN  lenses. 


OUTSIDE  CONE 


(a)  Terrell's  Lens  With  Conical  Outer  Surface 


(b)  Lens  Shapes  In  This  Thesis  With  Conical 
Inner  Surface. 


Figure  1.  Lens  Shapes 


II.  THE  HOMOGENEOUS  LENS 


A.  THEORY 

Snell's  law  is  the  cornerstone  of  contemporary  lens  de¬ 
sign  in  homogeneous  optical  materials.  In  the  HIN  case 
Snell's  law  is  used  in  the  familiar  form 

N^  sin  =  N2  sin  02  (1) 

where  N^  is  the  index  of  refraction  of  the  material  corres¬ 
ponding  to  surface  intercept  angle  0^,  with  respect  to  the 
surface  normal,  at  the  interface  between  surfaces.  Lens 
geometry  and  the  relative  values  of  Ni  determine  the  resultant 
optical  behavior. 

B.  ASSUMPTIONS  AND  SIGN  CONVENTION 

In  order  to  simplify  the  design  and  analysis  of  the 
seeker  lens  problem,  certain  assumptions  have  been  intro¬ 
duced.  Although  energy  loss  upon  transmission  through  the 
lens  at  each  surface  is  calculated,  it  is  assumed  that  the 
light  is  monochromatic  radiation,  time  dependent 
electric  and  magnetic  fields.  Light  impinging  upon  the  lens 
is  assumed  to  have  a  planar  wave  front  as  if  propagating  from 
an  object  at  infinity.  The  presence  of  a  shock  wave  attached 
to  the  lens  is  ignored  as  are  any  other  regions  of  expansion 
or  compression  in  the  flow  field  about  the  lens  [6] .  Further¬ 
more,  the  index  of  refraction  of  the  free  stream  is  assumed 
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to  be  equal  to  that  of  the  interior  space  behind  the  lens 
even  though  computer  routines  were  written  with  the  flexi¬ 
bility  to  process  unequal  values. 

The  sign  convention  used  is  a  right  handed  system  with 
spatial  coordinates  positive  to  the  right,  up,  and  out  of 
the  page  as  seen  by  the  reader.  All  angles  are  assumed 
positive  counterclockwise  from  point  of  reference  unless 
otherwise  noted.  All  linear  dimensions  are  implicitly  non- 
dimensionalized  with  respect  to  the  maximum  radius  of  the 
right  circular  cone  forming  the  inside  surface  of  the  lens. 

C.  HIN  LENS  DESIGN 

The  lens  design  procedure  consists  of  calculating  both 
the  loci  of  points  forming  the  outside  surface  and  the  slope 
of  the  surface  at  each  point  in  the  meridian  plane.  It  is 
convenient  to  approach  the  problem  by  placing  a  point  source 
of  light  at  the  design  focal  point  F  and  calculating  succes¬ 
sive  refracted  ray  paths  U'  using  Snell's  law  at  points 
on  the  inside  surface  in  Figure  2.  Points  E^  are  formed  by 
the  intercept  of  the  refracted  ray  and  the  slope  of  the  out¬ 
side  surface  as  extended  from  the  previous  point  E^_^.  The 
tangent  to  the  surface  which  will  refract  the  ray  in  the  de¬ 
sired  direction,  U,  may  again  be  found  from  Snell's  law. 

The  accuracy  of  the  calculated  surface  increases  as  the 
spacing  between  points  decreases  or  as  angle  AU"  becomes  very 
small.  The  final  ray  is  parallel  to  the  lens  axis  and  point 
D  corresponds  to  point  B,  which  is  the  origin.  The  area 
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formed  by  the  triangle  ABEn  is  opaque  to  prevent  rays  from 
scattering  through  the  opposite  side  of  the  lens. 


A  particular  lens  is  begun  by  specifying  the  focal 
length  BF,  the  cone  half-angle  a,  the  indices  of  refraction 
N^,  ,  and  N^,  the  cone  radius  at  the  edge  R,  the  ray 

direction  angle  U,  and  the  total  number  of  rays  to  be  traced. 
The  thickness,  T,  of  the  lens  at  the  edge  (E-^D^)  must  also 
be  specified  in  order  to  define  which  of  the  family  of  possi¬ 
ble  outer  surfaces  will  be  calculated.  In  general,  angle  U 
will  be  taken  to  be  zero  in  order  to  investigate  objects  on- 
axis  at  infinity.  R  will  always  be  set  equal  to  one  and  to 
nondimensionalize  all  linear  dimensions  every  length  is  im¬ 
plicitly  expressed  as  a  ratio  with  respect  to  R.  The  angle 
6  at  the  focal  point  is  measured  from  the  lens  axis  to  and 
is  expressed  as 


tam  0 


(BF  -  R* cotan  a] 


Thus,  if  1+1  is  the  total  number  of  rays  to  be  traced,  then 


S  -  J’AU' 


where  J  is  the  ray  number.  To  begin  tne  lens  design  the 
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coordiantes  x2(J),y2(j)  of  point  Dj  must  be  found.  Here 

the  subscript  2  refers  to  the  inside  surface  where  x^(j) 

y. (J)  are  the  coordinates  of  E_  on  the  outside  surface. 

-L  u 

Following  Kingslake  [7],  Q  in  Figure  3  may  be  expressed 
as 

Q  =  BF  sin  U"  =  x2(J)  sin  U"  +  y2(J)  cos  u"  (5) 


But , 


y2(J)  •  =  x2(j)  tan  a 


(6) 


so  that 

0 

BF  sin  U"  »  x2(J) [sin  U"  +  tan  a  cos  U"]  (7) 


Now,  x2(J)  and  y2(J)  are 

BF  sin  U" 

X2(J)  *  [sin  u" + tan  a  cos  U"] 


(8) 


(  -n  BF  in  U"  tan  a 

^2 1  *  fsin  U"+tan  cx  cos  U"] 


(9) 


Looking  now  to  find  the  incident  angle  which  the  ray  forms 
with  the  inside  surface  at  Dj  it  can  be  seen  that 

BF  -  x2(J) 

y2(J> 
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tan  (a  +  II) 


(10) 


Figure  3.  Ray  Trace  Geometry 


Substituting  into  Equation  (1),  Snell's  law,  the  refraction 


angle  I 2  can  now  be  written  as 

-1  N3 

I2  *  sin  A[(^)sin  l£]  (12) 

The  path  of  the  ray  inside  the  lens  may  now  be  expressed 
in  the  form  of  the  angle  U'  where 

U*  =  (J-a)  -  I2  (13) 

To  this  stage  all  rays  are  treated  the  same.  Now,  however, 
a  differentiation  must  be  made  between  the  outermost  ray 
which  defines  the  edge  of  the  lens  and  all  other  subsequent 
rays. 

For  the  outermost  ray  the  thickness,  T,  of  the  lens  at 
the  edge  can  be  specified  which  yields  x^l)  and  y^l) 
immediately : 


x-^  (  1) 

*  x2(l) 

-  T  cos 

U' 

(14) 

-  *2(1) 

+  T  cos 

U' 

(15) 

At  point  E^  we  may  write  Equation  (1)  again  to  obtain 

-1  ^1 

I,'  »  sin  ^[(rj^Osin  I.]  (16) 


where  angles  1^  and  I£  are  the  incident  and  refracted  angles 
at  the  outside  surface  and  are  not  yet  known.  Angles  1^ 
and  may  be  found  by  first  noting  that 


U  +  II  =  U'  + 


(17) 


then  substituting  Equation  (16)  into  (17)  and  rearranging 
terms  to  obtain 


-I  ^1 

*  U'  -  U  +  sin  [  (jp-)  sin  Ij 


(18) 


It  can  be  shown  after  some  algebraic  manipulation  that 
Equation  (18)  may  be  solved  for  the  angle  1^  in  the  form 


sin”1!- 


sin* (U'  -  U) 

T 


1/2 

-}  (19) 


[cos(U'-U)  -  N1/N2]  +  sin  (U1  -  U) 


Angle  1^  may  be  found,  if  desired,  by  substitution  of  Equa¬ 
tion  (19)  into  Equation  (16) .  Of  more  importance,  however, 
is  the  determination  of  the  slope  of  the  outer  surface  and 
the  surface  normal.  Now  that  angle  1^  is  known  these  slopes 
may  be  written  as 


=  SJLi 

dx  normal 


DYDXN(J) 


-  tan (U  +  Ix) 


(20) 


tti 


DYDXT(J) 


surface 

tangent 


cotan (U  +  L^) 


(21) 


Equations  (20)  and  (21)  may  be  evaluated  by  using  Equation 


(19)  for  I1. 


Now  that  all  the  parameters  are  known  for  the  first  ray. 


the  remainder  of  the  points  Ej  and  the  respective  slopes 


may  be  found.  Each  successive  ray  is  traced  as  before  by 
Equations  (4)  through  (13) .  Equations  (14)  and  (15) ,  how¬ 
ever,  may  not  now  be  used  since  the  lens  thickness  along  the 
ray  can  not  be  specified.  Instead,  the  intersection  of  the 
ray  and  the  surface  slope  from  the  previous  ray  is  used  to 
define  the  new  point  EJ+^.  The  intersection  is  found  by  first 
writing  the  equations  of  lines  representing  the  ray  and  the 
surface  tangent.  For  the  ray: 


-x  tan  U'  +  y2(J)  +  x.,(J)  tan(U') 


(22) 


For  the  surface  tangent: 


x  cot  (I^U)  +  y^J-1)  -  x^J-1)  cotd^U) 


(23) 


Equations  (22)  and  (23)  are  solved  simultaneously  to  yield 


the  coordinates  of  E_  which  are 

u 


y.(J)  - 


(24) 
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where : 


A  »  cot(i1+u)  cot(u') [y2(J)+x2(J) tan(U') ]  (25) 

B  »  y^J-1)  -  x^J-1)  cotdj^  +  U)  (26) 

C  -  1  +  cott^+U)  cot(U’)  (27) 

and 


xx(j)  »  cot(U')  t-y1(J)+y2(J)+x2(J)  tan(U')  1  (28) 


Now  Equations  (19) ,  (20) ,  and  (21)  are  used  to  calculate  1^ 
and  the  slopes  at  E_.  Therefore,  the  remainder  of  the  lens 
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surface  may  be  generated.  The  opaque  region  at  the  surface 
is  formed  by  extending  the  slope  of  the  surface  at  ER  to 
intersect  the  lens  axis.  Here,  K  refers  to  the  last  ray. 
The  nose  half-angle,  y,  thus  formed  is  given  by 


tan  y 


dx | surface,  ER 


(29) 


D.  SKEW  RAYS 

A  skew  ray  is  one  that  begins  from  an  off-axis  object 
point  and  enters  the  lens  either  in  front  of  or  behind  the 
meridian  plane  (z  *  0) .  For  every  skew  ray,  there  is  a 
corresponding  mirror  image  skew  ray  on  the  opposite  side  of 
the  meridian  plane  so  that  two  skew  rays  are  traced  at  the 


expense  of  only  one  calculation.  These  two  skew  rays  inter¬ 
sect  at  the  same  diapoint. 

Large  numbers  of  skew  rays  are  traced  through  a  lens 
in  order  to  study  lens  performance  at  different  obliquities. 
The  procedure  is  to  superimpose  a  grid  over  the  lens  aperture 
and  to  trace  rays  through  the  intersections  of  the  grid, 
through  the  lens  and  onto  the  image  plane.  An  image  plane 
spot  diagram  and  an  energy  density  plot  are  then  constructed 
for  study. 

To  accomplish  this,  the  aperture  grid  has  been  attached 
to  the  nose  (opaque  region)  of  the  seeker  lens  at  station  A; 
see  Figure  3.  The  plane  of  the  aperture  grid  is  tilted 
relative  to  the  lens  by  a  variable  angle.  A  transformation 
between  the  grid  coordinates  and  the  lens  coordinates  has 
been  derived  to  connect  skew  rays  from  grid  to  intercept  with 
the  outside  surface  of  the  lens.  Referring  to  Figure  4,  it 
can  easily  be  seen  that 

x'  =  [x  +  AB ]  cos  a  -  y  sin  a  (30) 

P  P 

y'  *  [x  +  AB]sin  a  +  y  cos  a  (31) 

P  P 

z  '  =  z  (32) 

where  is  the  tilt  of  the  grid  plane,  AB  is  the  length  of 
the  opaque  nose  portion  on  the  x-axis.  Both  z  and  z'  are 
positive  out  of  the  page. 


The  three  dimensional  lens  outer  surface  is  generated  by 
rotating  the  array  of  outer  surface  coordinates  in  the  merid 
ian  plane  through  2tt  about  the  lens  axis.  Each  pair  of 
coordinates  x^(J),y^(J)  thus  describe  a  circle  in  the  y-z 
plane  of  the  lens.  This  circle  transforms,  however,  to  an 
ellipse  in  the  grid  plane  given  by  the  equation 


y'  =  [x  +  AB]  sin  a  ±  VR2  -  z' 

D 


cos  a 


where  R  is  the  radius  of  the  circle  in  the  x-y  plane.  R 
may  be  expressed  by  the  familiar  equation  of  a  circle 


2  2 
=  y  +  z 


in  the  y-z  plane  of  the  lens,  or  by  solving  Equation  (33; 
for  R  in  the  grid  plane 

R2  =  z  ' 2  +  ^ - (x  +  AB)  tan  a  )  2  (35) 

cos  a_  p 

P 

With  the  aid  of  the  foregoing  groundwork,  the  x,y,z  coor¬ 
dinates  of  the  ray  intercept  with  the  outer  surface  may  be 
found.  To  see  how  this  is  accomplished,  first  refer  to  Figure 
5.  The  skew  ray  will  pass  outside  the  circle  formed  by  rotat¬ 
ing  some  point  x^(J),y^(J)  and  inside  the  next  circle  formed 
by  rotating  x^  ( J-l)  ,y1  ( J-l)  .  In  so  doing  the  skew  ray  will 
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Poin 


intercept  an  imaginary  cone  formed  by  the  two  circles. 
Now,  the  y',z'  of  the  skew  ray  in  the  grid  plane  and  a 


trial  value  (J)  from  the  array  of  surface  points  are  sub¬ 
stituted  into  Equation  (35)  and  a  value  for  R  is  calculated. 
If  R  is  greater  than  y^ ( J)  but  less  than  y^(J-l)  than  the 
appropriate  circles  have  been  found.  If  not,  another  trial 
<'alue  of  x1(j)  is  picked.  In  practice  the  x^(J)  corresponding 
to 


yx(J)  =  z'  (36) 

is  chosen  as  the  first  trial  value  and  subsequent  trial 
values  are  x^J-l),  x^J-2)  and  so  forth. 

The  equation  of  the  cone  passing  through  the  two  circles 
may  be  derived  by  again  referring  to  Figure  5  and  noting 
that 


QP  =  QP  ' 


n — i 

*  Vy  +  2  =  R 


(37a) 


QP  =  y  =  m[x  -  xx(J) ]  +  y1(J) 


(37b) 


where  the  slope  m  is  given  by 


^Recall  that  a  straight  line  between  points  Ej  was  used 
to  form  the  outer  surface  in  the  meridian  plane.  These 
approximations  are  valid  only  if  the  spacing  between  points 
or  circles  is  very  small. 


b 


y—mr  — —m  7^— t  1  ■  "J*  TJr 


p>: 


3  o? 


m  = 


y^J-1)  -  y^(J) 
x^J-1)  -  xx(J) 


(38) 


and  therefore 


7  y1  ( J-l)  -  y,  (J)  2 

R  *  tl4'c?=n--x^jrllx-xli>1  +  yi(J)) 


(39) 


By  equating  Equation  (39)  to  Equation  (35)  and  solving  for 
x,  the  expression  for  the  x-coordinate ,  x^ ,  of  the  outside 
surface  ray  intercept  may  be  found.  After  some  algebra,  this 
reduces  to  the  complicated  relationship 


PARI  +  ./PARI  PAR? 
PAR?  "  PAR?  *  PAR2 


where: 


PARI  *  ab  +  ce  -  b2AB  -  c2d 


2  2 

PAR2  =*  b  -  c 


and  where 


(40) 


(41) 


(42) 


PAR3  =  a2  +  b2  (AB)  2  -  2abAB  +  2  '  -  (e-cd)  2  (43) 


a  *  y'/cos  a 


(44) 


tan  a. 


(45) 
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y1(J-l)  -  y^J) 

c  =  Ix1(J-l)”3q‘(jJ1 

d  *  x^J) 

e  *  y^J) 


(46) 

(47) 

(48) 


It  has  been  found  that  the  plus  sign  in  Equation  (40)  gives 
the  correct  values  until  PAR1/PAR2  becomes  greater  than 
x^(l)  at  which  time  the  negative  sign  must  be  used.  It  now 
follows  from  Equations  (31),  (32)  and  '40)  that  the  yQ  and 
coordinates  of  the  intercept  are 


*0 


_ 

cos  ap 


(XQ  +  AB)  tan  otp 


(49) 


and 


Zq  =  z*  (50) 

The  next  step  in  tracing  the  skew  ray  is  to  ascertain  the 
direction  cosines  of  the  ray  inside  the  lens  after  refrac¬ 
tion  at  the  outside  surface.  This  is  accomplished  in  a 
series  of  steps  beginning  with  the  determination  of  the 
angle  <p  which  the  skew  ray  makes  with  the  normal  to  the 
surface  at  the  point  of  intercept.  Angle  <p  may  be  found  by 
taking  the  scalar  product  of  the  direction  cosines  of  the 
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skew  ray  and  the  surface  normal.  Thus  to  find  <p,  the  surface 
normal  must  first  be  obtained. 


As  might  be  surmised  by  the  reader,  the  direction  cosines 
of  the  surface  normal  can  be  found  by  taking  the  gradient 
of  Equation  (39)  with  Equation  (37)  substituted  for  R. 

Several  orders  of  magnitude  my  be  gained  in  accuracy,  how¬ 
ever,  if  the  normal  Nm  is  found  by  interpolation  in  the 
meridiem  plane  and  then  rotated  to  the  point  of  intercept  by 
angle  y  as  in  Figure  5.  This  interpolation  procedure  involves 
multiplying  the  ratio  S/ST  with  the  difference  between 
DYDXN(J)  and  DYDXN(J-l)  and  adding  the  product  to  DYDXN(J) . 

DYDXNpIM  =  ^(DYDXN(J-l)  -  DYDXN(J)  )  +  DYDXN(J)  (51) 


where  subscript  PIM  refers  to  point  of  intersection,  meridian. 
Values  of  the  slopes  of  the  normals  in  the  meridian  plane, 
DYDXN(J) ,  are  given  by  Equation  (20)  and  S  and  ST  are  the 
linear  separation  of  points  P",  PIM  and  P",  P  respectively. 
S/ST  is  given  by 


S 

ST 


2  2  1/2 

I  <R- yx(J) )  +  (xpiM“xi(J))  ] 

/A  2  .  2,1/2 

(Ax  +  Ay  )  ' 


(52) 


Therefore,  the  normal  vector,  Nm,  in  the  meridian  plane  is 


(DYDXN__„ 


)i  + 


j 


(53) 


and  since 


tan  y 


(54) 


the  normal  vector  at  the  point  of  intersection  on  the 
outside  surface  becomes 


N' 


1  *  -1  Z0  A  -1  zo  * 

(nvnxM - +  cos  (tan  —-)  j  +  sin(tan  ~)  k  (55) 

u u  PIM  y0  y0 


and 


o  1/2 

| N ' |  *  [DYDXNpIM”Z  +1]  (56) 

A 

Now  the  scalar  product  of  the  skew  ray  vector,  R,  and  surface 
normal,  N* ,  may  be  performed  to  find  $ .  Since  the  acute 
angle  between  these  vectors  is  required,  the  dot  product 
must  be  written 


R*N '  *  j  N '  |  cos U-<t>) 

Substituting  Equation  (55)  for  N'  and  noting  that 


Ki  +  Lj  +  Mk 


cos  a  i  -  sxn  a  j 
P  P 


the  left  side  of  Equation  (56)  may  be  expanded  to 


R*N ' 


008  an  -1  z0 

DYDXN  ‘  sin  ap  cos  (tan 

PIM  p  y0 


(57) 


(58) 


(59) 
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Solving  Equation  (56)  for  the  angle  $  and  introducing 
Equation  (58)  ,  the  expression  for  angle  <t>  becomes 


7T-COS 


-1 


COS  aD  -1  z0 

DYDXH^  -  “p  COS<tan 


-2  2  -1  Z0  -2  -1  Z0  1 
[DYDXN_tm  *  +  cos^  (tan  ^-)  +  sin  *(tan  i  —■)  ] 
PIM  Yn  Yn 


(60) 


Now  that  the  acute  angle,  <J>,  between  the  ray  and  the  outside 
surface  normal  is  known,  0',  the  angle  between  the  refracted 
ray  inside  the  lens  and  the  surface  normal  may  be  found  using 
Snell's  law. 


-1  1 1 

4>  *  ■  sin  (sr-  sin  <j>)  (61) 

w2 

The  next  step  is  to  ascertain  the  direction  cosines  K',  L', 

M'  of  the  refracted  ray  inside  the  lens.  Following  Kingslake 
[7] ,  Figure  6  shows  the  optical  vector  relationship  between 
R,  the  skew  ray,  R' ,  the  refracted  skew  ray,  N' ,  the  surface 
normal,  and  the  indices  of  refraction  of  the  two  media. 
Algebraically  this  relationship  is  written  as 

N2R'  *  N^R  +  (N2cos  r-N1  cos  4>)N'  (62) 

By  resolving  Equation  (61)  into  component  form,  the  direction 
cosines  of  R'  may  be  found.  After  rearranging,  K',  L',  and 


M'  become 


5T—  K  +  (  COS  $  1 
2 


rj-  COS  4>)k 
W2 


js—  L  +  (COS  <J>  ' 
w2 


rj—  COS  <J>)  l 
n2 


M'  =  —  M  +  (cos  <j> '  -  cos  4>)m 


where  K,  L,  and  M  are  given  by  Equation  (57)  and  k,  l,  m  are 


DYDXN, 


-1  Z0 
cos (tan  — ) 

y0 


-I  Z0 
sin  (tan  x 


Again,  |n' |  in  Equations  (66)  through  (68)  is  expressed  by 
Equation  (56)  . 

At  this  stage  in  the  process  of  tracing  the  skew  ray, 
the  coordinates  of  the  external  ray  intercept  xQ ,  yQ,  zQ 
are  known  as  are  the  direction  cosines  K',  L',  M'  of  the 
skew  ray  inside  the  lens.  The  next  step  is,  of  course,  to 
find  the  ray  intercept  with  the  inside  surface  (cone)  and 
the  direction  of  the  ray  subsequent  to  refraction.  The 
intermediate  steps  are  similar  to,  if  not  identical  with. 


the  foregoing.  One  major  difference,  however,  is  that  the 
inside  conical  surface  may  be  expressed  analytically  and  the 
intercept  coordinates  x^,  y^,  z^  may  be  found  exactly  without 
approximation.  Proceeding,  it  can  be  seen  that  if  x^,  y^, 
z ^  were  already  known,  Equations  (63) ,  (64) ,  and  (65)  could 
be  rewritten  as 


K 


.  _  xi  ~  x0 


D 

.  „  yi :  yo 


(69) 


(70) 


M1 


Zi  ‘  20 


(71) 


where  D  is  the  distance  between  surface  intercept  points. 
Upon  rearranging : 


xi  *  DK'  +  xQ  (72) 

y±  »  OL*  +  yQ  (73) 

Zi  ■  DM'  +  Zq  (74) 

Furthermore,  there  exists  a  relationship  between  x^,  y^,  and 
z ^  which  is  given  by  the  expression  for  the  inside  conical 
surface : 


z2  - 


2  2 
x .  tan  a 


0 


(75) 


where  a  is  the  half-angle  of  the  cone.  Substituting  Equa¬ 
tions  (72),  (73),  and  (74)  into  Equation  (75)  and  solving 
for  0,  it  is  seen  that 


D 


(76) 


where 

P1  3  L'y0  +  M'Z0  "  K'xQtan2a  (77) 

P2  =  L'2  +  M'2  -  K 1 2  tan2  a  (78) 


and 

2  2  2  2 
p3  “  Yg  +  Zq  -  xQ  tan  a 


(79) 


Correct  values  for  D  are  obtained  by  using  the  minus  sign  in 
Equation  (76).  Now  that  D  is  known,  values  for  K',  L',  M' , 

Xq ,  yQ ,  and  Zq  are  substituted  into  Equations  (72) ,  (73) ,  and 
(74)  to  yield  the  coordinates  of  the  inside  surface  intercept 
point. 

In  the  case  of  the  inside  surface,  the  gradient  may  be 
used  to  obtain  the  surface  normal  at  the  point  of  intercept 
since  this  surface  has  been  expressed  analytically.  If 
Equation  (75)  is  denoted  by  f(x,y,z)  then 


where  is  the  unit  vector  in  the  direction  of  the  surface 
normal  at  the  point  of  intercept.  Here, 


**  -  Hi* 


or  alternately. 


Vf  «  fxi  ♦  fyj  +  t%k 


(81) 


(82) 


Thus,  the  direction  cosines  of  the  normal  k',  l' ,  m'  may 
be  expressed  as 


k' 


V 


m' 


Vf 


Vf  I 

f 

_ z_ 

Vf  I 


(83) 


(84) 


(85) 


where : 


f 


2x-  tan^a 


(86) 


d 


V f  |  =  2  [x?  tan4  a  +  y?  +  z?] 


1/2 


(89) 


Now  that  the  two  unit  vectors  R*  and  N'  are  known  the  inci¬ 
dent  acute  angle  of  the  skew  ray  with  the  inside  surface, 

<{k  ,  may  be  ascertained  by  again  taking  the  scalar  product  as 
in  Equation  (57) .  Thus, 

4>i  =  TT-cos"1  [K'k'  +  L  +  M’m']  (90) 


and  using  Equation  (1)  again 


♦i 


sm 


-1 


N2 


(91) 


where  is  the  acute  angle  between  the  skew  ray  and  the 
surface  normal  after  refraction  at  the  inside  surface. 
Direction  cosines  K" ,  L",  M"  of  the  skew  ray  after  refraction 
are  found  analogously  to  K',  L',  and  M'  in  Equations  (63), 

(64)  and  (65).  Here,  however,  N2,  N3,  K',  L',  M',  4k,  4>^, 
and  k',  l*  ,  m'  are  substituted  for  ,  N2,  K,  L,  M,  $,  4> '  , 
and  k,  l,  m  respectively.  Thus 


K" 


K '  +  (cos  4>  • 
N  3  1 


N2 

rj—  cos  •  )  k  1 
N  2  i 


(92) 


L" 


n2  n2 

rj-  L'  +  (cos  4>  :  -  cos  . )  l ' 


(93) 


% 

>  . 

i 

f  * 

.  + 

i  m 

*  .  ’ .  . 


^  -w-  »-w  ■■■»  V.  "■  **-'  "- 


and 


M" 


N  N 

jj—  M'  +  (cos  ~  jj—  cos  4>j_)n»' 


(94) 


Finally,  the  skew  ray  intercept  with  the  image  plane  may 

be  found.  The  image  plane  is  treated  as  another  surface 

along  the  path  of  the  skew  ray,  and  the  intercept  coordinates 

x.  ,  y .  ,  z.  are  easy  to  find.  Equations  (72),  (73)  and 
im  J  lm  im 

(74)  may  be  used  again  in  the  form 


x 


im 


P'K"  +  xA 


(95) 


im 


*  D'L"  +  y. 


Z 

im 


D'M"  + 


(96) 


(97) 


where  D',  in  this  instance,  is  the  linear  separation  between 
coordinates  x^,  yi,  z^  at  the  inside  surface  and  x^m  ,  Yim/ 
z of  the  image  plane.  Furthermore, 


x 


im 


BF 


(98) 


is  the  equation  of  the  image  plane.  It  follows  that 

(BF  -  x.) 

D’  =  - TTii  ~  (99) 
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after  substitution  of  Equation  (98)  into  Equation  (95) . 
Therefore, 


(100) 


(101) 


Hence,  the  skew  ray  has  been  traced  onto  the  image  plane 

at  the  focal  point  of  the  lens.  The  y.  coordinate  of  the 

corresponding  mirror  image  skew  ray  is  the  same  as  Equation 

(100).  The  z.  coordinate,  however,  is  the  negative  of 
un 

Equation  (101)  since  the  mirror  image  ray  is  behind  (when 
viewed  along  the  z-axis)  the  meridian  plane.  After  tracing 
a  complete  set  of  skew  rays  through  the  lens  a  spot  diagram 
may  be  plotted.  Clearly,  the  number  of  rays  to  be  traced 
depends  entirely  upon  the  incremental  size  of  the  aperture 
grid  chosen.  According  to  Kingslake  [7] ,  at  least  100  rays 
must  be  traced  to  give  a  fair  approximation  of  the  actual 
image.  In  addition  to  the  spot  diagram,  an  energy  density 
plot  may  now  be  constructed  by  counting  the  number  of  rays 
within  progressively  larger  radii  from  the  image  centroid  and 
then  plotting  the  number  of  rays  as  a  function  of  radius. 
Here,  each  ray  is  assumed  to  contain  a  unit,  nondimensional , 
an  amount  of  radiant  energy  for  convenience. 


E.  RADIANT  ENERGY  LOSS 

In  reality,  each  ray  loses  intensity  upon  transmission 
at  each  interface.  Of  the  total  amount  of  energy  contained 
in  each  ray,  a  fraction  IT  will  be  transmitted,  a  fraction  IR 
will  be  reflected,  and  a  fraction  IA  will  be  absorbed  by  the 
medium  into  which  the  ray  is  propagating.  Since  it  has  been 
assumed  that  absorption  is  negligible,  it  must  be  true  that 

IT  +  IR  =  1  (102) 

Furthermore,  the  relative  amounts  of  transmitted  and  reflected 
electromagnetic  energy  may  be  calculated  by  the  well-known 
Fresnel  Equations  which  state  the  dependency  of  IT  and  IR 
upon  the  angle  of  incidence  and  the  indices  of  refraction 
at  the  interface.  That  IT  and  IR  are  further  dependent  upon 
the  orientation  of  the  electric  vector  with  respect  to  the 
geometry  of  ray  incidence  is  fundamental  to  the  boundary 
conditions  which  govern  the  form  of  the  Fresnel  relations  as 
derived  in  Hecht-Zajac  [8] .  Since  this  thesis  examines  lens 
response  to  monochromatic  radiation,  the  indices  of  refraction 
are  not  considered  as  a  function  of  wavelength;  further,  the 
electric  vector  orientation  is  assumed  to  be  rapidly  and 
randomly  changing  with  time.  By  time  averaging  field  com¬ 
ponents,  it  may  be  seen  that  the  reflectance  is 

IR  =  |(r  l  +  rjj)  (103) 
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where 


rX 


cos  $  -  (N^ 
- ji 

COS  $  +  (Nfc^ 


2A. 

-  sm  $) 

- ~2 — 

-  sm  $) 


1/2 

T7I 


(104) 


and 


M2  .  ,„2  .  2a  .1/2 

Nfci  cos  <J>  -  (Nti  -  sm  <Jk) 

~2 - 2 - 5 — T7T~ 

Nti  cos  ^  +  ^Nti  ”  s^n  ^  ' 


(105) 


Here,  is  the  ratio  of  the  index  of  refraction  of  the 

transmission  side  of  the  interface  to  the  index  of  refrac¬ 
tion  of  the  incident  side.  From  Equation  (102)  it  now  follows 
that  the  transmittance  through  the  interface  is 


(106) 


The  total  transmittance  through  the  lens  is  simply  the  product 
of  Ifc  at  the  outside  surface  with  that  of  the  inside  surface 
where  4k  is  substituted  for  4>  in  Equations  (104)  and  (10  5)  . 

Total  internal  reflection  of  the  ray  may  occur  at  the 
inside  surface  if  the  incident  angle  becomes  too  large. 
Following  Reference  (7) ,  this  occurs  when  angle  4k  is  equal 
to  or  greater  than  tt/2.  Thus,  Snell's  law  becomes 

N3 

sin  4>.  *  —  (107) 

i  «2 

and  any  ray  with  <Jk  equal  to  or  greater  than  this  will  be 
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totally  internally  reflected.  In  this  thesis  such  rays  are 
labeled  "failed  rays"  since  they  fail  to  intersect  the  image 
plane.  If  N2  is  1.5  and  is  1.0,  the  incident  angle  for 
total  internal  reflection  is  41.81°,  or  greater. 

F.  OPTICAL  PATH  LENGTH  (OPL) 

The  optical  path  length  of  a  skew  ray  is  an  analytical 
tool  with  which  the  researcher  may  ascertain  the  phase  of  a 
ray  at  the  end  of  the  path.  By  so  doing,  the  image  diffrac¬ 
tion  pattern  may  be  constructed  which  shows  the  addition  or 
subtraction  of  amplitude  depending  upon  relative  phase.  Since 
each  ray  must  begin  with  the  same  phase,  monochromatic  radia¬ 
tion  is  used  for  diffraction  experiments.  OPL  is  included 
here  only  as  a  matter  of  interest.  The  calculation  of  optical 
path  length  is  simply  the  svun  of  the  geometrical  path  segments 
of  a  ray  multiplied  by  the  corresponding  index  of  refraction 
of  the  medium  for  that  segment. 


III.  THE  GRIN  LENS 


A.  THEORY 

This  thesis  assumes  a  spherically  symmetric,  inhomogene¬ 
ous,  isentropic  medium  in  which  the  refractive  index  varies 
from  point  to  point  but  is  independent  of  direction  at  each 
point.  The  refractive  index  is  a  function  of  the  coordinates 
of  the  points  of  the  region  being  considered.  The  problem 
of  describing  the  resulting  curved  paths  of  rays  in  such 
GRIN  materials  has  been  solved  long  ago  in  the  form  of  a 
single  second  order  vectorial  differential  equation.  Marchand 
[9]  has  shown  that  the  solution  to  the  differential  equation 
in  the  case  of  spherical  gradients  can  be  written  in  polar 
coordinates  in  the  plane  of  the  ray  as 


9  =  9„  +  e 


77r 


dr 

1 


(108) 


Here  rQ  and  9g  are  values  of  r  and  9  at  a  convenient  reference 
point  on  the  ray  and  |e|  is  a  scalar  constant  along  the  ray 
given  by 


e  =  e  n0r0sin*0 


(109) 


Referring  to  Figure  7,  r  is  measured  from  the  center  of 
symmetry  of  the  index  function;  angles  9  and  0Q  are  measured 


counterclockwise  from  the  lens  axis;  angles  and  are 
measured  counterclockwise  from  the  radius  vector  to  the 
instantaneous  ray  direction  vector,  e  is  the  sign  function 
expressed  as 


e  =  sgn  U0|  *  ±1 


(110) 


where 


l^o  I  <  tt/2  (111) 

In  Equation  (109)  ,  n^  refers  to  the  index  function  n  evaluated 
at  rQ .  The  form  of  n  has  been  chosen  as 

n  *  n(0  .r)  (112) 

s 

where  O  denotes  the  position  of  the  center  of  symmetry  on 
s 

the  lens  axis  (shown  in  Figure  7)  and  r  is  the  radial  coor¬ 
dinate  from  Os  to  the  point  in  question.  This  form  allows 
the  study  of  the  effects  on  lens  performance  as  Os  is  changed. 

More  specifically,  in  order  to  allow  the  analytical 
integration  of  Equation  (10  8)  ,  the  expression  used  for 
n(Og,r)  is  a  generalized  version  of  that  employed  by  Luneburg 
as  described  by  Marchand  [9] .  Here 

r  r  «  2 , 1/2 

n  «  [a  +  b (-k—  )  ] 


(113) 


whereas  Luneburg  used  specific  values  for  a  and  b.  Generaliz¬ 
ing  the  index  function  enables  the  strength,  or  percent  change 
of  the  gradient  to  be  varied  as  well  as  the  algebraic  sign 
of  the  gradient.  If  the  parameter  b  is  negative,  a  decreasing 
parabolic  gradient  results;  conversely  a  positive  b  yields 
an  increasing  parabolic  gradient.  By  visualizing  a  plane 
wave  front  passing  through  a  spherical  gradient,  it  may  be 
seen  that  the  negative  gradient  results  in  light  rays  bending 
forward  at  the  center,  as  shown  in  Figure  7,  whereas  a  posi¬ 
tive  gradient  has  the  opposite  effect.  This  principle  was 
used  by  Wood,  as  related  by  Marchand  19],  in  constructing 
simple  lenses  having  plane  faces  and  a  radial  index.  The 
Wood  lens  acted  as  a  converging  or  diverging  lens  depending 
upon  the  sign  of  the  gradient  used. 

Additionally,  it  must  be  noted  that  if  the  parameter  b 
is  equal  to  zero.  Equation  (113)  reduces  to 

n  -  /a  (114) 

which  is  a  HIN  lens  having  constant  index  of  refraction. 

This  fact  has  facilitated  the  correlation  of  GRIN  and  HIN 
computer  trace  algorithms. 

Equation  (108)  may  be  integrated  by  a  change  of  variable 


using  the  relation 


|  sin”1 

2e2/r2  -  a 

-  sin  1 

r  2e  /ro " a  i) 

/  9  “5  0 

t  J—  ”2  '  2 

'  w  Va*  +  4be*/R:  J 

L  Va z  +  4be  VR;J  ’ 

Equation  (116)  gives  0  as  a  function  of  r.  This  equation 
can  be  easily  solved  for  r  as  a  function  of  6  in  the  form 


ja+Wa  tip-Sin|-2£(9-e0)tsin’1|-— :  .  ■  ■■  , 

(  V  RZ  L  Va2+4bei/R 


in"1  - 


■Jr 


(117) 

The  instantaneous  direction  of  the  ray  at  any  point  r,9 
may  be  ascertained  by  using  the  invariance  of  e.  Hence, 


e  =  enr  siniji  *  enn  rn  sin 


(118) 


1  n0r0sin  * 
sm  [ - — - J 


(119) 


Furthermore,  the  orientation  of  the  plane  of  the  ray  may  be 
easily  deduced  since  every  ray  in  a  spherical  medium  is  a 
plane  curve  lying  in  a  plane  through  the  center  of  symmetry 


5 


Using  this  fact,  Marchand  [9]  has  shown  that  a  suitable 
conversion  from  coordinates  r  and  9  in  the  plane  of  the  ray 
to  global  Cartesian  coordinates  may  be  written  in  the  form 


X 

s 

r  (6 

X0 

^  +  n  po> 

(120) 

y 

= 

r  (6 

p-  +  n  q0) 
r0  u 

(121) 

z 

r  ( 6 

+  n  a  ) 

(122) 

ro  0 

Here  pQ,  qQ,  JlQ 

are 

the 

initial  direction 

cosines  of  the 

ray  at  r0,90;  xQ 

'  *0 

|f  and  Zq  are  the  Cartesian  coordinates 

corresponding  to 

V 

• 

O 

CD 

The  parameters  6 

and  n  are  given 

by  Marchand  as 

n 

= 

sin 

9/sin  <pQ 

(123) 

6 

a 

cos 

9  -  n  cos  ipQ 

(124) 

It  should  be  noted  that  Equation  (117)  may  become  singu¬ 
lar  for  certain  rays  where  <|/  or  ipg  become  very  close  to  zero 
or  tv.  This  singularity  may  be  more  easily  seen  in  Equation 
(123)  where  n  becomes  indeterminate  as  both  9  and  ipQ  approach 
zero  and/or  it.  In  practice  these  conditions  occur  when  0 

3 

is  located  either  far  out  in  object  space,  coincident  with 
B,  or  on  the  image  side  of  the  lens;  see  Figure  8.  Positioning 
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the  center  of  symmetry  at  these  locations  has  therefore 
been  avoided. 

B.  ASSUMPTIONS 

The  aforementioned  assumptions  for  the  HIN  case  also 
apply  here.  Additionally ,  it  must  be  assumed,  as  depicted 
in  Figure  7,  that  the  GRIN  lens  could  or  soon  may  be  fabri¬ 
cated  from  a  sphere  of  dielectric  material  with  the  required 
spherically  symmetric  parabolic  gradient. 

C.  GRIN  LENS  DESIGN  PARAMETERS 

In  the  HIN  lens,  the  available  design  parameters  are 
basic.  These  include:  F,  R,  T,  a,  U,  and  N2-  Parameters 
available  for  varying  the  design  of  the  GRIN  lens,  however, 
include  those  of  the  HIN  case  but  expand  the  index  of  refrae 
tion  variable  N2  into  0g,  a,  +b,  and  -b.  These  additional 
lens  design  parameters  greatly  expand  the  lens  designer's 
power  to  bend  radiant  energy  to  his  will. 

D.  GRIN  LENS  DESIGN 

The  GRIN  lens  design  procedure,  although  paralleling 
that  of  the  HIN  case,  is  somewhat  more  complicated  in  that 
the  rays  are  now  curved  and  che  index  of  refraction  varies. 
Accordingly,  the  same  design  process  is  used  but  with  more 
intermediate  calculations  required. 

The  additional  calculations  arise  since  the  intercept 
of  the  GRIN  ray  with  the  surface  tangent  cannot  be  solved 
in  closed  form,  and  an  iterative  solution  must  be  used. 


c 


The  reader  will  note  that  Equation  (117)  will  yield  r  if 
(9-9q)  is  known.  Hence  the  iterative  procedure  is  to 
"guess"  (9-9q)  based  on  the  HIN  coordinates  which 

are  calculated  as  in  the  homogeneous  lens.  Thus 


9h  "  8t  8o 


(125) 


V* 


in  Figure  9  is  used  in  Equation  (117)  to  find  point  p.  The 
tangent  to  the  ray  path  at  point  p  is  now  extended  to  inter¬ 
cept  the  surface  tangent  again  using  the  homogeneous  inter¬ 
cept  relations  to  find  The  prime  superscripts 

indicate  successive  iteration  values.  The  distance  d  is 
employed  as  a  measure  of  the  error  of  point  p.  If  d  is  not 
within  an  acceptable  margin  then  9^  is  calculated  based  on 
xlH,ylH  and  the  Procedure  repeated  to  find  p'  in  Figure  9. 

If  d'  is  not  within  allowable  error  then  the  iteration  con¬ 
tinues  until  it  is  acceptable.  In  practice,  this  iteration 
procedure  has  proved  to  be  extremely  rapid,  rarely  requiring 
more  than  three  iterative  steps  before  converging.  Slight 
modifications,  however,  must  be  introduced  to  handle  a  posi¬ 
tive  gradient.  Furthermore,  if  0g  is  located  outside  the 
lens  proper,  a  decreasing  angular  increment  must  be  subtracted 
from  9h  to  ensure  that  each  radial  vector  rH,  r^,  ...  inter¬ 
sects  the  ray  during  iteration.  These  modifications  to  the 
intercept  iteration  procedure  are  recorded  in  the  program 
listing  for  program  GISL  (for  Gradient  Index  Seeker  Lens) . 

Refer  to  Appendices  A,  B,  and  C  for  a  full  description  of  GISL. 
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Figure  9.  Grin  Ray  Intercept  Geometry 


With  the  overall  iterative  intercept  procedure  now  being 
clear,  the  mathematical  details  of  the  GRIN  lens  design 
follow.  Referring  again  to  Figure  9,  the  GRIN  ray  is  traced 
as  in  the  HIN  design  procedure  to  obtain  the  coordinates 
x2(J),y2(J)  on  the  inside  surface  in  the  meridian  plane.  Now, 
before  Snell's  law  can  be  used,  the  index  of  refraction  must 
be  found  at  these  coordinates.  Using  Equation  (113) 


rA  ,  1/2 


n2  =  [a  +  b(^)2] 


(125) 


where 


rQ  =  \/ (x2(J)  +  OB)  2  +  y2(J) 


(126) 


and  where 


J  (x  "(1)  +  OB)  2  +  y,(l)  2 


Rz  =  v  (x2(l)  +  OB) ~  +  y2 


(127) 


Furthermore,  by  inspection  it  can  be  easily  seen  that 


*o  =  "  -  U'  -  90 


(128) 


U'  is  identical  to  the  HIN  case  and 


'Here  Equation  (125)  must  be  substituted  for  N2  in 

.it.on  (12) . 


(129) 


9o  tan  [y2(j)  +  obj 

In  Equations  (126)  through  (129)  OB  is  the  line  segment  from 

0  to  B  and  is  defined  as  a  positive  quantity  to  the  left 
s 

of  B  for  algebraic  ease  of  manipulation.  By  substitution 
into  Equation  (109)  the  scalar  constant,  e,  of  the  ray  can 
now  be  found.  Furthermore,  since  angle  U'  is  now  known 
coordinates  x^'y^H  ntay  calculated  using  Equations  (14) 
and  (15)  or  Equations  (24)  through  (28) . 

For  the  first  ray  defining  the  edge  of  the  lens.  Equa¬ 
tions  (14)  and  (15)  are  employed,  and,  since  no  iteration 
is  required,  angle  (0,j,-0q)  is  used  in  Equation  (117)  to  find 
r  immediately.  Here 


cos  0T 


X1H  +  0B 


'H 


(130) 


and 


■H 


n/Tx 


1H 


2  2 
+  OB )  +  y 


1H 


(131) 


Thus 


x.  ( 1)  =  r  cos  9_  -  OB 


(132) 


yxu) 


r  sin  8, 


(133) 


are  the  coordinates  of  the  first  point  on  the  outside 
surface.  Since  e  is  known.  Equation  (119)  is  employed  to 
find  the  angle  i|>.  Angle  c  in  Figure  10  is  employed  to 
translate  ip  into  the  ray  direction  with  respect  to  the  lens 
axis  by  the  relation 


?  =  it  -  (<l<  +  9t> 


(134) 


Therefore,  by  substituting  £  for  U*  in  Equations  (19) ,  (20)  , 
and  (21),  DYDXN(l)  and  DYDXT(l)  may  be  found.  Here  again,  N2 
in  Equation  (19)  must  be  replaced  by  n2  as  given  by  Equation 
(113)  evaluated  at  x^d)  ,y^(l)  . 

Now  that  the  first  point  on  the  outside  surface  is  known, 
along  with  the  surface  tangent,  the  remainder  of  the  K  number 
of  rays  may  be  processed  to  yield  the  balance  of  the  outside 
surface.  Each  subsequent  intercept  between  ray  and  tangent 
must  be  iterated.  Thus,  unlike  the  first  ray,  once  x^H,  y  , 
rH'  rp'  an<*  ^  are  known,  x1{J  1  and  y1H  ’  are  ascertained  by 
the  substitution  of  the  coordinates  of  point  p  for  x2(J),y2(J) 
in  Equations  (24)  through  (28)  with  U’  replaced  by  £.  The 
coordinates  of  point  p  are 


Xp  =  rp  COS  eT  "  0B 


v  =  r  sin  9_ 


(135) 


(136) 


Therefore 

w 

=  {A  +  B} 

1  c  1 

where 

A  » 

cot(I^  +  U)  cot 

e  [yp  +  xp  tan  z  ] 

B  = 

yx ( J-l)  -  x1(J 

-1)  cot  [I1  4-  u] 

C  = 

1  4-  cot  (I1  4-  u) 

cot  z 

and 

T  _  =5 

sin2  (C-U) 

LcosU-U)  -N1/n2J  +  sin  (C-U) 


Note  that  n2  in  Equation  (140)  is  found  from  Equation  (113) 

evaluated  at  r  . 

P 

Additionally, 


X1H '  *  cot  +  yp  +  xptan;] 


(1411 


Now  the  error,  d,  may  be  evaluated  as 


Here  d  is  compared  to  1  *  10  .  If  d  is  larger  than  this 
value  then  the  entire  procedure  is  repeated  by  substituting 

x1h'  'ylH '  ^or  xlH,ylH  an<*  so  ^ort^*  Once  the  error  criteria 
are  satisfied 

x,  (J)  =  x  (143) 

1  P 

tw«  •  • 

yL(J)  *  yp  (144) 

and  the  next  ray  is  processed.  Correlation  between  GRIN 
(with  b  set  equal  to  zero)  and  HIN  design  procedures  run 
with  identical  parameters  has  shown  agreement  to  the  fifth 
and  sixth  decimal  places. 

E.  SKEW  RAYS  IN  GRIN 

GRIN  skew  ray«  are  handled  analogously  to  the  homogeneous 

case  with  the  same  coordinate  transformation  from  grid  plane 

to  global  coordinates  being  required.  It  is  only  after  the 

initial  directions  cosines  of  the  GRIN  skew  ray  K‘,  L',  M' 

are  found  that  the  differences  between  GRIN  and  HIN  appear. 

The  only  exception  to  this  being  the  use  of  Equation  (113) 

in  Snell's  law  for  refraction  at  the  interface.  Since  the 

GRIN  skew  rays  display  curvature  in  a  plane  through  0  ,  K', 

s 

L',  and  M'  are  constantly  changing  until  intercept  with  the 
inside  conical  surface.  Therefore,  not  only  must  the  plane 
of  the  skew  ray  be  analytically  described,  but  the  final 
values  of  K',  L',  M'  must  be  found.  Due  to  the  nature  of 


GRIN  rays,  the  procedure  for  finding  the  ray  intercept  with 
the  conical  inside  surface  is  different  from  both  the  HIN 


case  and  the  iteration  procedure  employed  in  the  meridian 
plane  due  to  the  multiplicity  of  the  geometry  encountered. 
The  Newton- Raph son  iteration  routine  has  been  found  to  be 
ideal  for  this  purpose. 

To  begin,  the  magnitude  and  direction  of  the  initial 
radius  vector  r^  from  0g  for  Xq ,  yQ,  zQ  must  be  ascertained. 
The  magnitude  is  given  by 

1/2 

r„  =  [<x„  +  OB)2  +  y*  +  ,145) 

Therefore  the  unit  vector  in  the  direction  of  the  intercept 
is 


A 


*0  *  0B 


/V 

i 


(146) 


The  plane  of  the  ray  may  be  fully  described  by  the  vector 

/S  A 

normal  to  the  plane.  Two  vectors,  r q  and  R  lie  in  the  plane 

A 

of  the  ray.  R  is  the  unit  vector  in  the  initial  direction 
of  the  ray  after  refraction  and  described  using  direction 
cosines  as 


R  =  K '  i  +  L’j  +  M'k  (147) 


Thus,  the  plane  of  the  ray  may  be  described  by  the  cross 
product 


[I 


rQx  R 


(148) 


or 


N 


N  i+N  j+N  k 
p0x  p0y  p0z 


where 


N  =  rn  M'  -  rn  L* 

pox  °y  02 


N  =  rrt  K'  -  r„  M' 
PQy  Oz  Ox 


and 


N 


p0: 


r0xL'  "  r0yK ' 


(149) 


(150) 


(151) 


(152) 


In  Equations  (150),  (151),  and  (152)  rQx,  rQy,  and  rQz 

/v 

refer  to  the  x,  y,  and  z  components  of  rQ  in  Equation  (146) . 

/x  /\ 

Furthermore,  the  angle  ip ^  between  Tq  and  R  may  be  found  from 
the  dot  product  as 


♦ 0  =  cos"ltr0xK'  +  r0yL’  +  r0zM'>  (153) 

Now  the  scalar  invariant  e  may  be  found.  Substituting  known 
values  into  Equation  (109) : 


e  =  e  [a  +  b(^-)  “]  rQ  sin 
z 


where  Rz  is  unchanged  from  that  found  during  previous  calcu¬ 
lations  for  the  lens  shape  by  Equation  (127)  . 

With  the  foregoing  groundwork  established,  the  intercept 
of  the  ray  with  the  inside  surface  may  be  calculated.  The 
Newton-Raphson  iteration  scheme  requires  the  calculation  of 
the  radius  vector  from  0g  to  the  cone  by  geometrical  methods 
and  the  radius  to  the  ray  by  GRIN  theory.  The  difference 
between  the  two  radii  is  then  divided  by  the  difference 
between  the  derivatives  of  the  two  functions.  The  resulting 

quantity  is  subtracted  from  the  trial  angle,  9  ,  in  the  plane 

P 

of  the  ray,  to  give  a  new  trial  angle  0^.  The  process  is 
continued  until  the  difference  between  radii  is  less  than 
lx  10~5.  The  first  trial  angle  is  measured  to  a  reference 
HIN  intercept  as  if  the  material  were  homogeneous  since  the 
actual  GRIN  ray  curves  only  slightly.  The  coordinates  of 
this  HIN  intercept  point  are  designated  x^H,  y^H,  ziH  and 
are  derived  using  the  HIN  equations  as  before.  To  obtain 
the  first  trial  angle,  the  scalar  product  between  rQ  and 

A 

r.„  is  used.  Here 

In 


2  2  2 

=  [  (x,„  +  OB)  +  y , „  +  z  .rrl 


(155) 


(156) 


n  - 


riH 


XiH  *  0B  :  +  ^iH  ^  +  fiH  J 
riH  riH  riH 


Thus, 


cos"1[rOxriHx  +  rOyriHy  +  cOzrLHz] 


(157) 


where  r^Hx,  riHy'  an<*  riHz  are  x'  y'  an<*  z  comPonents 
of  respectively. 

In  Equation  (117)  8^  is  substituted  for  (9-0q)  to  yield 
r  as  required  by  the  iteration  procedure.  The  geometrical 
radius,  rg,  is  not  so  easily  acquired.  First,  note  that  the 
equation  of  the  plane  of  the  ray  inside  the  lens  is  given  by 


(  e? 


N 


P0X(X-V  +  NF0y(^V  +  Np02(2‘Z0)  * 


(158) 


i 
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Secondly,  the  equation  of  the  conical  surface  is  given 
by  Equation  (75) .  The  combination  of  the  plane  of  the  ray 
and  the  cone  yield  the  loci  of  possible  intercept  points  on 
the  inside  surface.  In  Cartesian  coordinates,  the  sum  of 
Equations  (75)  and  (158)  is 


x(Np0x'xtan  ]  +  yCNpoy  +  y]  +  z[Npoz  +  zl  -  Npoxxo 


Npoyyo  Npozzo  ® 


(159) 


Equation  (159)  must  be  transformed  into  coordinates  r  and  0 


90 


B 


1  «? 


in  the  plane  of  the  ray.  The  transformation  is  made  possible 
using  Equations  (120)  through  (124)  .  Upon  substitution, 
and  after  solving  for  r^,  Equation  (159)  becomes 


where 


-B2 


,  V4T 


4A~  A2 


2  2  2  2 
A2  =  Bi  +  c±  ~  A^  tan  a 


+  N  z. 
poz  0 


(160) 


(161) 


B2  =  BlNpoy  +  ClNpoz  +  AlNpox  +  2Al0Btan"a  (162) 


C_  =  OBN  +  OB2tan2a  +  N  x.  +  N  y_ 
2  pox  pox  0  poy-'O 


(163) 


1 


and 


~-(xn  +  OB)  +  nK* 
rQ  u 


(164) 


B. 


SYr 


+  nL' 


(165) 


6  z, 


+  nM' 


(166) 


Here,  5  and  n  are  found  from 


n  =  sin  9  /sin  ipQ 
P  ^ 


(167) 
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h 


S  m 


and 


<5  =  cos  9  -  cos  ipn 

p  0 


(168) 


which  follow  from  Equations  (123)  and  (124)  .  The  plus  sign 

in  Equation  (160)  yields  the  correct  values.  Now  that  r  of 

the  ray  and  of  the  surface  are  known,  the  derivatives  of 

r  and  r  with  respect  to  9  at  9  must  be  found.  It  can  be 
3  P 

shown  that  for  the  ray 


dr 

d9 


where 


— — ^  cos  [-2  e  9  +  sin  ^  (A-)  ] 

P  3 


2e 


(169) 


2e2/r2  -  a 
^a2  +  4be2/R2 


(170) 


The  derivative  of  rg  is  somewhat  more  complicated.  With 
persistence,  however,  it  can  be  shown  that 


dr 

~d9 


i  dB-  dA  ,  B2  C  '1/2 

771  tA2  ~de~  "  b2  'dF]  +  2  t??!  +  aT1  tc3] 
2A2  4A^  2 


(171) 


where  A2,  B2,  and  C2  are  given  by  Equations  (161)  through 
(168)  and 


1  ?  dB_  _  dA- 

l^(A2  b2  inr  "  b2  a2  -dT) 


C2  ^2. 

I1  ^ 

A2 


(172) 
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(173) 


Therefore,  the  r  and  0^  to  the  intercept  are  found  by 
iteration  of  revised  trial  values 


In  practice  the  quotient  of  differences  in  Equation  (173) 
is  reduced  by  a  factor  of  1.3  to  slow  convergence  and  provide 
stability.  The  number  of  iterations,  however,  rarely  exceeds 
five . 

Values  for  r  and  0  are  now  transformed  into  Cartesian 

P 

coordinates  x^,  y^,  by  substitution  into  Equations  (120) 
through  (124) .  Next,  the  values  of  the  direction  cosines 
K',  L',  and  M'  at  the  intercept  are  needed.  The  angle  \p 
between  the  radius  vector  and  the  tangent  to  the  ray  at 
intercept  may  first  be  deduced  from  the  scalar  invariant, 
e .  Hence 

<j>  =  sin  1(— | — )  (174) 

n2rPI 

Where  e  is  known,  is  evaluated  at  rpi;  rpi  is  the 
radius  to  the  intercept  as  found  by  the  iteration  above. 

Three  constraints  on  the  direction  cosines  may  be  written. 
These  are: 

1)  The  scalar  product  of  the  radius  vector,  rpi,  and  the 

A 

instantaneous  ray  direction  vector,  R,  at  intercept. 

2)  The  scalar  product  of  the  normal  to  the  plane  of  the 

/\  A 

ray  N^q  and  R. 
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3)  The  sum  of  the  squares  of  the  direction  cosines  must 
sum  to  unity. 

Mathematically,  the  above  constraints  are  written  as 


(x.  +  OB)  y.  z. 

- = -  K’  +  — L'  +  — = -  M' 

T  r  t 

pi  pi  pi 


=  cos 


N„_  K'  +NrtL'+N.M'  =  0 

POx  pOy  pOz 


K 


2  2  2 
^  +  L ' £  +  M’  = 


(175) 


(176) 


(177) 


To  find  K’,  L',  and  M',  Equations  (175),  (176)  and  (177)  are 
solved  simultaneously.  It  can  be  shown  that  the  solution 
1 eads  to 


M'  = 


-B6 


b; 


4A 


2  A, 


(178) 


L'  =  A5  -  B5M' 


K’ 


where 


2  2 
(1  -  L  '  *  -  M '  ; 


1/2 


(179) 


(180) 


A,  =  U  *  ( 


^POz  2  2  _  2  Np0y  Np0z  B 

N  n  1  5  v2  °5 

P°x  NP0x 


N 


+  (Jj^-)  2  +  1 

Nd0x 


(181) 
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2  N  -  N  -  N  _  - 

B  =  _ P°Y  P°2  A  -  1  [1  +  (.,.£P-2)  21A  B 

6  N  n  A5  J  u  n  ;  A5“5 

pOx  pOx 


(182) 


[l  +  (-E°2-)2]a^  -  1 

N  «  5 

pOx 


(183) 


NpOxcos* 


N  _  B.  -  N  n  A. 
pOx  4  pOy  4 


(184) 


N  n  C .  -  N  _  A. 
pOx  4  pOz  4 


N  .  B.  -  N  .  A. 
pOx  4  pOy  4 


(185) 


Furthermore,  in  Equations  (184)  and  (185),  A4 ,  B4 ,  and  C4 
are  the  coefficients  of  K',  L'  and  M'  in  Equation  (175)  . 

Now  that  the  direction  cosines  of  the  ray  are  known  at 
the  point  of  intersection  with  the  inside  surface,  the  index 
of  refraction,  n^,  is  computed  by  substituting  rpI  into 
Equation  (113) .  The  remainder  of  the  skew  ray  trace  to  the 
image  plane  is  identical  to  the  homogeneous  procedure. 
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and 
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-  Np0xC4  ~  Np02A4 
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Furthermore,  in  Equations  (184)  and  (185),  A4,  B^,  and 
are  the  coefficients  of  K',  L'  and  M*  in  Equation  (175)  . 

Now  that  the  direction  cosines  of  the  ray  are  known  at 
the  point  of  intersection  with  the  inside  surface,  the  index 
of  refraction,  n2,  is  computed  by  substituting  rpI  into 
Equation  (113) .  The  remainder  of  the  skew  ray  trace  to  the 
image  plane  is  identical  to  the  homogeneous  procedure. 


LENS  PERFORMANCE  PARAMETERS 


The  function  of  the  seeker  lens  is  to  focus  electromag¬ 
netic  energy  either  reflected  from  or  emitted  by  the  target 
onto  a  detector.  Angular  displacement  of  the  target  with 
respect  to  the  missile  body  axes  as  well  as  target  angular 
rate  information  are  both  desired  outputs  from  the  seeker. 
Hence,  it  is  not  only  important  just  to  be  able  to  detect 
the  target  by  focusing  energy  into  a  spot  on  the  detector, 
it  is  equally  important  that  this  spot  be  as  small  as  possi¬ 
ble  to  enable  the  precise  position  of  the  spot  on  the  detec¬ 
tor  to  be  discerned. 

The  ability  of  a  lens  to  focus  an  object  to  a  small  spot 
does  not  guarantee  the  quality  of  the  image.  For  a  FLIR 
optical  system  extensive  effort  is  expended  to  obtain  an 
image  with  minimum  aberration  [10].  Seeker  optics,  however, 
are  generally  non-imaging  devices  where  the  pressure  of  the 
different  aberrations  does  not  detract  from  the  function  of 
the  seeker  as  long  as  a  tight  image  is  maintained  [11] . 

Accordingly,  the  most  important  parameter  by  which  seeker 
lens  performance  is  judged  is  that  of  spot  size  at  different 
obliquities.  Since  the  image  found  at  the  focal  point  is  not 
necessarily  circular  nor  equally  dense,  the  standard  devia¬ 
tion  in  the  y  and  z  directions  with  respect  to  image  centroid 
is  used  to  define  spot  size. 

Therefore,  once  the  coordinate  pairs  of  all  the  skew  rays 


have  been  calculated  in  the  image  plane,  the  first  step  in 


the  analysis  of  lens  performance  is  to  calculate  the  image 

centroid  [12].  If  the  spot  diagram  is  composed  of  N  rays, 

the  centroid  location,  y  is 

c 


(186) 


Since  there  is  symmetry  about  the  x-y  plane,  the  z  coor¬ 
dinate  of  the  centroid  will  always  be  zero. 

The  next  step  in  finding  the  spot  size  is  to  find  the 
standard  deviations  of  the  spot  diagram  in  the  y  and  z 
directions.  This  is  accomplished  by  summing  the  squares  of 
the  differences  of  the  intercept  coordinates  with  respect 
to  the  centroid  and  then  dividing  by  the  number  of  rays. 

The  standard  deviations  then,  are  given  by 


(187) 


(188) 


(189) 
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Of  further  interest  in  appraising  lens  performance  is 
the  energy  density  of  the  image  as  a  function  of  radius  from 


the  centroid.  In  nondimens ional  form, this  is  simply  the 
number  of  rays  in  the  spot  diagram  within  a  succession  of 
circles  of  increasing  size  overlaid  about  the  centroid.  Here 
each  ray  is  assumed  to  carry  a  unit  amount  of  radiant  energy. 
This  type  of  plot  facilitates  the  comparison  of  different 
lens  designs  by  detailing  the  distribution  of  energy  within 
each  image.  Clearly,  it  is  desirable  to  have  as  much  energy 
as  possible  concentrated  very  close  to  the  image  centroid. 
Between  two  lenses  with  equal  spot  sizes,  the  preferred  lens 
has  more  energy  concentrated  within  a  smaller  radius. 

Every  ray,  however,  does  not  deliver  an  equal  amount  of 
energy  to  the  focal  plane.  It  is  prudent,  therefore,  to 
include  as  a  performance  parameter  the  average  ray  intensity. 
Again,  for  N  rays,  the  intensity  I  of  each  ray  is  summed  and 
normalized  by  N  to  yield 
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Finally,  each  lens  design  is  checked  for  "failed  rays”. 
The  reader  will  recall  that  these  rays  fail  to  intercept  the 
image  plane  due  to  total  internal  reflection,  total  external 
reflection,  or  failure  to  intercept  the  inside  surface  within 
the  bounds  of  the  lens.  Hence,  a  lens  design  with  fewer 
"failed  rays”  or  no  "failed  rays"  at  all  is  a  preferred  lens. 


Although  there  are  many  other  performance  criteria  by 
which  lenses  are  compared,  the  foregoing  parameters  are  more 
than  sufficient  to  judge  the  merit  of  preliminary  seeker 
lens  designs.  It  should  be  noted,  however,  that  notwithstand 
ing  the  fact  that  the  image  centroid  and  standard  deviations 
were  used  as  stepping  stones  to  obtain  image  spot  size,  they 
have  significant  meaning  of  their  own.  The  standard  devia¬ 
tions  o  and  <?z  inform  the  lens  designer  as  to  the  horizontal 
and  vertical  spread  of  the  image.  Image  centroid  location, 
y  ,  at  increasing  obliquities  is  of  obvious  importance  since 
excessive  displacement  will  cause  the  image  to  miss  the 
detector  entirely  and  would  dictate  the  necessity  for  a 
second  lens  element  to  dampen  the  movement.  Furthermore, 
since  Line  of  Sight  (LOS)  measurement  accuracy  to  the  target 
is  highly  dependent  upon  the  linearity  of  y  as  a  function  of 
the  lens  tilt  angle,  lens  designs  which  exhibit  a  greater 
degree  of  such  linear  behavior  are  the  preferred  designs. 

To  summarize,  it  is  sufficient  to  note  that  although 
spot  size  is  the  most  important  of  the  performance  parameters 
every  other  parameter  has  a  significant  impact  on  the 
performance  of  a  particular  lens. 


RESULTS  FOR  THE  HOMOGENEOUS  LENS 


The  performance  of  the  homogeneous  lens  is  presented 
primarily  as  a  comparison  with  which  to  compare  the  performance 
of  the  GRIN  lens.  Here,  the  relationship  of  spot  size  to 
increasing  lens  obliquity  and  lens  thickness  are  presented 
as  well  as  the  image  centroid  movement  as  a  function  of 
obliquity.  Additionally,  the  reader  is  introduced  to  the 
four  basic  computer  plots  used  to  display  the  results:  lens 
shape,  object  plane  with  superimposed  skew  ray  grid,  image 
plane  or  spot  diagram,  and  energy  density.  These  plots  were 
generated  on  the  VERSATEC  Plotter  using  arrays  of  data  points 
produced  by  program  GISL  on  the  IBM  3033  mainframe  computer. 

To  begin.  Figure  11  shows  the  homogeneous  lens  shape. 

The  first  of  the  four  basic  plots,  the  lens  shape  plot,  pre¬ 
sents  the  lens  side  view  in  the  meridian  plane.  The  outer 
surface  (curved)  and  the  inner  surface  (cone)  are  constructed 
by  connecting  the  points  x^(j)  ,  y1(H)  and  x2(J) ,y2(J) > 
respectively,  by  straight  lines.  At  the  apex  of  the  lens 
about  the  lens  axis  is  a  trapezoidal  region  which  represents 
the  opaque  nose  area.  In  the  legend  are  listed  the  lens 
design  parameters  and  the  significant  calculated  dimensions 
of  the  lens.  Since  all  linear  dimensions  are  implicitly 
normalized  with  respect  to  the  maximum  inside  radius,  R, 
the  lens  may  be  scaled  up  or  down  by  multiplying  each  dimen¬ 
sion  by  a  factor  of  Rnew/R*  Beginning  at  the  top,  the 
parameters  listed  in  the  legend  are  explained  in  Table  1. 
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TABLE  1 


PARAMETER 

TYPE  OF  PARAMETER 

MEANING 

F/# 

DESIGN 

F  number.  F/#  *  F/2R 

F 

design 

* 

Focal  length  from  B. 

OSYMB 

DESIGN 

OB.  Line  segment  Q  to  B 
(positive  to  left) . 

Immaterial  in  HIN 

R 

DESIGN 

* 

Maximum  radius  of  cone 

ALPHA 

DESIGN 

* 

a--cone  half-angle,  radians 

T 

DESIGN 

* 

Edge  thickness 

U 

DESIGN 

Incident  ray  offset  angle 
(design) *,  radians 

I 

DESIGN 

Number  of  iterations. 

1+1  *  number  of  rays 

N1 

DESIGN 

Free  stream  index  of 
refraction 

N3 

DESIGN 

Index  of  refraction  of 
interior  lens  cavity 

A 

DESIGN 

a  in  n-(r)  *  >/a  +  b(r^/R^, 
gradient  refractive 
index  function.  N~  =  /a 
in  HIN. 

B 

DESIGN 

b  in  n2(r) .  Zero  in  HIN 

beff 

CALCULATED 

b  effective  *  b/R^ 

%  Index 

DESIGN 

Percent  change  in  n_(r) 
from  r  *  0  to  r  *  Rg 

%  Inside 

CALCULATED 

Percent  change  in  n2(r) 
along  inside  surface  from 
lens  axis  to  edge 

%  Outside 

CALCULATED 

Percent  change  in  n2(r) 
along  outside  surface 
from  opaque  region  to  edge 
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Table  1  (Continued) 


PARAMETER  TYPE  OF  PARAMETER  MEANING* 


%  Across 

CALCULATED 

Percent  change  in  n«(r) 
across  lens  from  lens 
axis  to  outside  surface 
at  the  thickest  point 

STATION  A 

CALCULATED 

x-cogrdinate  at  nose  of 
lens* 

GAMMA 

CALCULATED 

Y— nose  half-angle  gf 
opaque  nose  region, 
radians 

TSURFL 

CALCULATED 

Total  outside  surface 
length  from  Station  A 
to  the  edge 

* 

Refer 

to  Figure  2  for 

clarification 

The  HIM  lens  shape  has  a  convex  outer  surface  with  maxi¬ 
mum  thickness  on  axis  of  almost  ten  tiroes  the  edge  thickness. 
Although  the  lens  has  a  good  aerodynamic  shape  resembling  an 
ogive,  the  outer  surface  is  not  a  circular  arc  nor  can  a 
single  analytical  function  be  fitted  to  the  array  of  points 
describing  the  surface.  Note  that  the  nose  half-angle,  y, 
is  almost  identical  to  the  cone  half-angle,  a. 

All  lenses  have  been  designed  with  a  cone  half-angle  of 
45°  which  is  approximately  the  maximum  angle  for  which  aero- 
dynamically  efficient  lens  shapes  may  be  designed,  consider¬ 
ing  a  free  stream  Mach  Number  not  to  exceed  three.  Without 
exception,  overall  lens  performance  is  more  severely  degraded 
as  angle  a  is  reduced. 

Table  2  explains  the  legend  of  Figure  12,  which  is  the 
second  basic  plot.  Here  the  lens  is  depicted  as  seen  from 
the  skew  ray  grid  plane.  The  lens  tilt  angle,  a  ,  causes 

r 

the  equally  spaced  (in  J)  circles  descriging  the  surface  of  the 
lens  to  appear  as  ellipses.  In  Figure  12,  the  grid  spacing  has 
been  reduced  from  0.1,  which  is  normally  used,  to  0.3  to 
allow  identification  of  individual  rays  for  correlation  with 
the  image  plane  spot  diagram.  Although  the  small  number  of 
rays  used  is  not  sufficient  to  given  an  accurate  definition 
of  spot  size,  the  number  is  sufficient  to  describe  where 
rays  in  the  object  plane  are  being  focused  in  the  image  plane 
by  the  HIN  lens.  The  skew  rays  in  Figure  12  have  been  numbered 
in  the  order  in  which  they  wer^  proces  *d.  Actually,  only 
rays  1  through  19  were  actually  *  .*ac«d;  20  through  33  are 


TABLE  2 


Explanation  of  OBJECT  PLANE  Plot  Legend 


PARAMETER  TYPE  OF  PARAMETER 


MEANING 


ALFAP  Analysis  otp.  Lens  tilt  angle, 

radians 

Rays  Analysis  Total  number  of  rays 

processed 

Failed  Rays  Analysis  Total  number  of  rays 

failing  to  pass  through 
the  lens* 


* 

Failed  Rays  are  indicated  on  the  plot  by  a  diamond 
superimposed  on  the  grid  location  of  the  ray. 

mirror  image  skew  rays.  Ray  30  corresponds  to  ray  14,  for 
example.  Here,  the  lens  has  been  tilted  by  0.4  radians,  or 
22.9  degrees,  and  33  rays  have  been  processed  of  which  none 
have  faV ed  to  intercept  the  image  plane.  The  staircase 
pattern  has  been  added  in  this  case  in  order  to  show  the 
resulting  distortion  present  in  the  image  plane  (Figure  13) . 

The  Spot  Diagram  in  Figure  13  is  an  example  of  the  third 
basic  computer  plot;  see.  Table  3.  Unlike  most  Spot  Diagrams, 
this  example  has  the  individual  rays  numbered  for  comparison  with 
Figure  12;  also  the  resulting  distorted  staircase  pattern 
is  sketched.  By  cross-referencing  individual  rays  between 
Figures  12  and  13,  it  is  possible  to  recognize  where  certain 
areas  of  the  lens  are  focusing  rays  in  the  image  plane.  Rays 
1,  4,  20,  and  22  about  the  opaque  nose  region  form  a  coma 
tail  which  contributes  most  of  the  image  spread.  Rays  33, 
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TABLE  3 


ilanation  of  Spot  Diagram  Legend 


PARAMETER 

TYPE  OF  PARAMETER 

MEANING 

ALFAP 

Analysis 

oip,  lens  tilt  angle, 
radians 

SPOT  SIZE 

Performance 

or,  see  Equation  (189) 

Y-CENTROID 

Performance 

yc,  y-coordinate  of 
image  centroid;  see 
Equation  (186) 

SIGMA  Y 

Performance 

Oy ,  y-standard  deviation; 
see  Equation  (188) 

SIGMA  Z 

Performance 

a2,  z -standard  deviation; 
see  Equation  (187) 

RAY  COUNT 

Performance 

Number  of  rays  striking 
the  image  plane 

18,  and  19  from  the  bottom  portion  of  the  homogeneous  lens 
are  imaged  at  the  top  and  are  widely  separated  from  the  core 
of  the  image.  In  general,  the  regions  of  the  lens  which 
have  been  found  to  contribute  the  bulk  of  the  widely  spaced 
rays  are  the  immediate  nose  region  and  the  lower  portion  of 
the  lens . 

Furthermore ,  the  distortion  present  in  the  image  of  the 
staircase  pattern  clearly  shows  that  regions  closest  to  the 
nose  yield  the  greatest  distortion.  Horizontal  lines  are 
switched  end  for  end  and  tilted  approximately  45  degrees. 

The  upper  portion  of  the  lens  performs  the  best.  Rays 
3,  7,  11,  and  15  in  the  meridian  plane  are  focused  on  the 
yIM  axis  in  a  fairly  tight  region;  rays  10,  14,  17,  27,  30, 
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and  32  about  the  upper  periphery  in  the  object  plane  are  all 
focused  within  the  image  core. 

Table  3'  explains  the  legend  of  the  spot  diagram.  It  is 
seen  that  the  standard  deviation  in  the  y-direction  is 
approximately  five  times  that  of  the  z-direction.  This 
elongation  of  the  image  is  not  readily  evident  in  the  spot 
diagram  since  the  ordinate  and  abscissa  have  not  been  plotted 
with  equal  increments.  This  results  from  the  great  disparity 
between  Spot  Diagrams  of  the  various  lens  designs  studied. 

It  is  important,  therefore,  for  the  reader  to  take  careful 
note  of  the  relative  sizes  of  the  y_„  and  z_„  axes. 

Of  primary  importance  is  the  spot  size  in  Figure  13.  The 
value  for  spot  size  is  adversely  affected  by  the  poorest 
performing  regions  of  the  lens.  Were  it  not  for  these 
errant  rays,  the  spot  size  would  be  considerably  smaller. 

The  image  intensity  pattern  is  benefitted,  however,  by  the 
fact  that  the  rays  spread  the  farthest  from  the  centroid 
contribute  significantly  less  energy  per  ray  than  those 
being  focused  in  the  core  of  the  image.  Table  4  lists  the 
relative  intensities  of  the  primary  skew  rays  plotted  in 
Figure  13.  It  is  seen  that  the  high  intercept  angles  experi¬ 
enced  by  the  rays  closest  to  the  bottom  of  the  lens  and,  to 
a  lesser  extent,  those  near  the  nose,  result  in  higher 
reflectivity  and  lower  transmission  through  the  lens. 

In  the  legend  of  Figure  14 ,  the  average  of  ray  intensity 
is  given.  The  Spot  Diagram  Energy  Density  distribution  may 
be  seen  at  a  glance.  The  fraction  of  energy  (number  of  rays) 
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a)  Nuraberical  Order 


Ray 

Intensity 

1 

0.746 

2 

0.789 

3 

0.921 

4 

0.892 

5 

0.877 

6 

0.811 

7 

0.920 

8 

0.917 

9 

0.904 

10 

0.870 

In  order 

of  Descending 

Intensity  Ray(s) 

0.921 

3 

0.920 

7 

0.919 

11 

0.918 

12,15 

0.917 

8 

0.916 

16 

0.910 

13,17 

0.904 

9 

Ray  Intensity 


11 

0.919 

12 

0.918 

13 

0.910 

14 

0.887 

15 

0.918 

16 

0.916 

17 

0.910 

18 

0.640 

19 

0.576 
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0.640 
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0.576 
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shown  as  a  function  of  radius  from  the  Spot  Diagram  centroid 
is  normalized  with  respect  to  the  total  number  of  rays 
traces  (Rays,  in  Figure  12)  regardless  of  whether  all  of 
the  rays  successfully  intercepted  the  image  plane.  Thus 
the  fourth  basic  plot  may  be  used  in  conjunction  with  the 
Spot  Diagram  to  further  define  the  image  concentration  with 
respect  to  the  centroid. 

The  response  of  the  HIN  lens  to  increasing  tilt  angles 
is  given  in  Figures  d-1  through  D-2  3  where  the  index  of  refrac¬ 
tion  of  the  lens  has  been  set  at  1.5.  Whereas  the  spot  size 

at  a  *  0.0  is  very  small,  that  at  a  =  0.7  radians  is  quite 

P  P 

large  at  67%  of  the  lens  radius.  Figures  15  and  16  summarize 

spot  size  and  centroid  locations  for  the  lens.  Spot  size 

growth,  although  somewhat  irregular  at  the  higher  angles, 

is  pronounced.  Furthermore,  the  centroid  movement  is  seen 

to  be  an  approximately  linear  function  of  ctp  until 

a  =  0.3  and  easily  exceeds  the  radial  dimension  of  the  lens 
P 

at  higher  tilt  angles. 

Figure  17  shows  that  the  HIN  lens  may  be  slightly  improved 
by  increasing  the  edge  thickness,  T,  by  a  small  amount.  The 
lens  may  be  otherwise  tuned  to  improve  performance  at  certain 
tilt  angles  by  designing  the  lens  with  U  slightly  greater 
than  zero.  These  performance  improvements  are  practically 
insignificant,  however,  and  neither  lens  tuning  by  the  param¬ 
eter  U  nor  T  produce  improvement  across  the  spectrum  of  tilt 
angles.  Instead,  an  improvement  at  one  ap  usually  has  resulted 
in  a  degradation  at  others. 
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Figure  16.  HIN  Lens  Centroid  Movement  versus 
a  at  N.  *  1.5 


Finally,  the  homogeneous  lens  is  shown  with  index  of 
refraction  of  three  at  the  intermediate  tilt  angle  of  0.3 
radians  in  Figures  D-24  through  D-27.  The  lens  shape  required 
to  accommodate  the  higher  index  of  refraction  is  seen  to  be 
thinner  and  displays  less  outside  surface  curvature  than 
the  HIN  lens  with  Nj  =  1.5.  Spot  size  is  significantly  re¬ 
duced.  Since  infrared  lenses  such  as  germanium  generally 
have  fairly  high  refractive  indices,  improved  lens  per¬ 
formance  at  these  higher  values  is  encouraging. 


The  performance  of  the  GRIN  lens  is  similar  in  many 
respects  to  the  homogeneous  lens.  The  relationship  between 
characteristic  regions  of  the  GRIN  lens,  such  as  the  lower 
lens  portion,  and  where  these  regions  image  bundles  of  rays 
is  identical  to  the  HIN  lens  as  depicted  in  Figures  11  and 
12.  The  growth  of  spot  size  with  increasing  ap  and  the 
respective  movement  of  image  centroid  typical  of  the  homogene¬ 
ous  lens  is  clearly  displayed  by  the  GRIN  lens  as  well. 

The  measure  of  the  superiority  of  the  GRIN  lens,  there¬ 
fore,  lies  in  the  successful  correction  or  improvement  of 
the  deficiencies  seen  in  the  HIN  lens.  Here,  the  reduction 
of  spot  size  is  of  primary  concern. 

As  a  modest  example  of  the  ability  of  the  gradient  refrac¬ 
tive  index  to  reduce  spot  size,  the  GRIN  lens  design  shown 
in  Figure  18  is  examined.  This  lens  is  very  similar  in  shape 
to  the  homogeneous  lens  with  Nj  *  1.5.  Note,  however,  that 
unlike  the  HIN  lens,  the  nose  half-angle,  y,  is  slightly 
larger  than  that  of  the  cone,  indicating  more  outside  surface 
curvature.  The  object  plane  for  this  lens,  at  a  *  0.4,  is 

It 

shown  in  Figure  19.  Here,  the  error  in  focusing  rays  present 
in  the  image  plane  has  been  superimposed  over  the  grid  plane. 
yT1_  coordinate  error  contours  (with  respect  to  the  centroid) 
are  shown  in  the  left  half  of  the  plane  and  zIM  contours  in 
the  right.  These  contours  vividly  show  that  the  largest 
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errors  in  both  y  and  z  stem  from  the  nose  and  bottom  of  the 
lens.  The  best  performance  is  contributed  by  the  upper 
central  region.  The  reader  will  note  that  negative  errors 
in  z  occur  below  the  y^  =  0  plane  of  the  lens  and  that  nega¬ 
tive  errors  in  y  straddle  the  yQ  =  0  plane. 

The  corresponding  intensity  contours  are  displayed  in 
Figure  20.  Of  significance  here  is  that  regions  which  perform 
relatively  well  in  y  and  z  error  also  perform  well  in  trans¬ 
mitting  energy.  Regions  of  the  lens  which  are  characterized 
by  relatively  high  angles  of  incidence/  therefore,  perform 
the  poorest  and  have  the  most  to  gain  from  better  combinations 
of  gradient  index  and  Og 

The  Spot  Diagram,  Figure  21,  displays  an  image  pattern 
virtually  identical  with  the  HIN  lens  in  Figure  11.  The 
GRIN  spot  size,  or,  however,  is  slightly  smaller  by  approxi¬ 
mately  2%  and  both  a  and  a  are  correspondingly  smaller 

y  z 

which  indicates  superior  performance  by  the  GRIN  version. 

Image  centroid  for  the  GRIN  is  displaced  further  than  that 
for  HIN.  The  average  intensity  transmitted  by  the  GRIN  lens 
compares  favorably  at  0.88  as  opposed  to  the  average  inten¬ 
sity  of  the  HIN  at  0.85.  Further  evidence  of  superiority 
of  the  GRIN  lens  is  seen  in  the  form  of  a  steeper  slope  to 
the  Encircled  Energy  Plot  (Figure  22)  than  that  of  the  HIN 
lens  in  Figure  13. 

Although  the  example  GRIN  design  displays  marginal 
superiority,  some  GRIN  designs  do  not.  In  particular,  it 
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Figure  20  .  Intensity  Contours  of  the  Example  GRIN  Lens  Shown  in 


-IMRGE  PLRNE 


SPOT  DIAGRAM 


fiLFAP  *  0.400000 
SPOT  SIZE  =0.4785753 
T  CENTRQy)  =  -1.5782780 

5IGMP  T  «  0.  1868387 
S I GMR  Z  =  0.0421957 
RflY  COUNT  =33 

*  * 


Figure  21.  Example  GRIN  Spot  Diagram  at  otp  *  0.4 
Radians  for  GRIN  Lens  Design  Shown  in 
Figure  18 
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has  been  found  that  as  the  center  of  symmetry  of  the  index 
function,  0_,  is  moved  further  away  from  the  interior  of  the 
GRIN  lens,  performance  deteriorates. 

Two  series  of  GRIN  lens  designs  have  been  examined  at 
the  intermediate  tilt  angle  of  0.3  radians  in  order  to  define 
the  spectrum  of  lens  performance.  Both  design  series  explore 
the  GRIN  lens  for  0g  inside,  outside,  and  far  outside  the 
interior  of  the  lens.  The  first  series.  Figures  E-l  through 
E-92,  are  GRIN  designs  in  the  "low  range"  of  refractive  index 
with  the  index  parameter  "a"  set  at  2.25.  This  first  set  is 
compared  to  the  HIN  lens  with  N^  *  /a  *  1.5.  The  second, 
or  "high  range"  series  is  for  a  =  9.0,  as  compared  to  the 
HIN  at  N2  *  3.0,  and  are  displayed  in  Figures  F-l  to  F-128. 
Both  the  high  and  low  range  lenses  exploit  gradient  changes 
of  5,  10,  25,  and  50  percent  both  positive  and  negative, 
where  possible.  A  negative  50%  gradient  in  the  low  range 
is,  of  course,  not  possible  since  an  index  of  below  1.0 
would  result.  The  center  of  symmetry  was  not  located  at 
x  ^  0  due  to  the  aforementioned  singularity  encountered  at 
very  small  angles. 

The  lens  shapes  of  the  GRIN  designs  shown  vary  widely. 

The  outside  surface  may  be  either  convex  or  concave  and  ex¬ 
hibits  a  thinner  profile  at  higher  values  of  refractive 
index.  Although  all  of  the  resultant  GRIN  lens  shapes  are 
superior  aerodynamically  to  the  hemispherical  seeker  lenses 
currently  in  use,  the  convex  version  of  the  lens  has  more 
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obvious  applications  for  a  Sidewinder  type  tactical  missile, 
whereas  the  concave  lens  shape  is  more  applicable  to  the 
diffuser  of  a  ramjet  with  nose  inlet. 

Object  plane  diagrams  show  that  for  lens  designs  exhibit¬ 
ing  smaller  spot  size,  fewer  or  none  of  the  skew  rays  failed. 
This  fact  underscores  the  success  of  GRIN  in  controlling  image 
deterioration  contributed  by  the  nose  region  of  the  lenses. 
Although  the  lower  lens  and  nose  regions  still  create  the 
largest  of  the  image  spread,  the  extent  is  reduced. 

Figures  23  and  24  are  contour  diagrams  summarizing 
spot  size  performance  of  GRIN  lens  designs  in  the  low  and 
high  range  respectively.  From  these  diagrams,  it  can  be  seen 
immediately  that  the  best  performance  is  obtained  from  lenses 
having  the  center  of  symmetry  inside  the  lens.  From  this 
fact,  it  may  be  deduced  that  large  changes  in  refractive 
index  are  desired  along  the  surface  of  the  lens  rather  than 
across.  Except  for  isolated  regions  at  -5%,  increasing 
gradients  produce  smaller  spot  sizes.  The  positive  gradient 
in  the  high  range,  however,  has  an  almost  constant  spot  size 
from  +5%  to  +50%  where  an  improvement  of  only  0.6%  is  seen. 

In  both  the  low  and  high  ranges  the  positive  gradient  at 
50%  and  at  OB  =  0.05  has  proved  to  be  the  best  performing 
combination.  Again,  as  in  the  HIN  lens  the  high  range 
exhibits  the  smaller  spot  size. 

The  best  GRIN  lens  is  now  further  examined  over  the  range 
of  obliquities  up  to  0.8  radians.  Furthermore,  since  it  has 
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Figure  23.  Contour  Plot  Summary  of  GRIN  Lens  Spot  Size  Performance 
for  a  -  0.3  radians,  a  =  2.25 
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Figure  24 .  Contour  Plot  Summarizing  GRIN  Lens  Spot  Size  Performance 
at  a_  =  0.3  radians,  a  =  9.0 


been  determined  that  placing  the  focal  point  at  x  »  2.0, 
improves  spot  size  performance  at  lower  values  of  ap,  another 
series  of  plots  is  given  for  F/l.  With  the  focal  point  at 
x  *  1.5,  severe  degradation  results;  for  F  ■  3.0  only  modest 
differences  with  F/2  are  evident.  Appendices  G  and  H  show 
the  performance  of  the  "best"  GRIN  lens  at  F/l  and  F/2. 

Here  it  may  be  fully  recognized  that  the  resulting  be¬ 
havior  of  a  "good"  gradient  refractive  index  seeker  lens  is 
far  superior  to  that  of  the  HIN  lens  for  op  <_  0.6  radians; 
see  Table  5.  The  reader  will  note  that  there  are  no  failed 
rays  even  at  0.8  radians  (45.8  degrees).  Although  centroid 
movement  at  F/2  is  slightly  greater  than  the  HIN  lens,  by 
changing  to  F/l,  this  can  be  corrected;  see  Figure  25.  The 
F/l  version,  moreover,  displays  a  range  of  spot  sizes  signi¬ 
ficantly  smaller  than  that  for  the  F/2  lens  below  ap  ■  0.54 
radians  or  the  HIN  lens  below  0.68  radians;  see  Figure  26. 

That  the  lens  displays  more  desirable  performance  with  a 
shorter  focal  length  is  a  surprising  but  highly  desired  re¬ 
sult.  Since  seeker  optics  systems  are  volumn  limited,  any 
reduction  in  the  lens  focal  length  is  beneficial  to  the  final 
design  and  packaging  requirements . 

It  is  significant  for  the  reader  to  note  that  centroid 
movement  beyond  a  value  of  ±1.0  exceeds  the  physical  dimen¬ 
sions  of  the  lens  radius.  If  the  interior  of  the  lens  mount 
is  of  similar  dimension,  it  follows  that  the  requirement  for 
^CENTR  to  ke  ^ess  than  *1.0  restricts  the  F/l  lens  to  approxi¬ 
mately  0.65  radians  of  tilt  whereas  the  F/2  lens  is  further 
restricted  to  less  than  0 . 3  radians .  The  superior  performance 
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Figure  25.  "Best"  GRIN  Lens  Centroid  Moment  for 
F/Numbers  of  F/l  and  F/2 


SPOT  SIZE  VS.  LENS  TILT  ANGLE 


TABLE  5 


Comparison  of  HIN  and  GRIN  Lens  Designs 


Performance  Parameter  or  Feature 

HIN 

GRIN* 

Shape  of  Outer  Surface 

Convex 

Concave 

Fraction  of  Failed  Rays 

at  a  »  0.7  radians 

P 

0.14  3 

0.0 

Linearity  of  Centroid 

Motion;  a_  for  Ay/y 
deviation*of  10% 

0 . 4  radians 

0 . 8  radians 

Spot  Size  for  ap 

a  *0.1  radians 

D 

0.190 

0.029 

a  *0.2  radians 

P 

a  *  0.3  radians 

D 

0.247 

0.442 

0.063 

0.112 

mi 

a  *  0.4  radians 

□ 

0.466 

0.201 

Cr 

a,  *  0.5  radians 

D 

0.637 

0.313 

a  *  0.6  radians 

P 

0.632 

0.490 

a  *0.7  radians 

D 

0.667 

0.688 

ap  *  0.8  radians 

N/A 

0.947 

Spot  radius  for  80% 
of  energy  §  a  *  0.4 
radians  ^ 

Does  not 
attain  80% 
within  0.50 

0.28 

Average  relative  intensity 
of  skew  rays  @  a  *  0.4 
radians  ** 

* 

0.82 

0.44 

* 

"Best"  GRIN  lens  with  50%  positive  gradient,  0  *  0.05, 

in  the  F/l  configuration.  s 
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of  the  F/2  lens  above  0.54  radians  is,  therefore,  unusable. 
Attempts  to  improve  upon  this  performance  by  increasing  the 
edge  thickness,  T,  or  by  slightly  adjusting  the  parameter  U, 
as  with  the  homogeneous  lens  only  produced  uniformly  degraded 
performance  in  every  respect. 

Although  the  "best"  GRIN  lens  resulted  from  a  positive 
50%  gradient  change  in  refractive  index,  it  has  been  noted 
that  this  configuration  was  only  slightly  better  than  the 
same  lens  with  a  positive  5%  change.  Furthermore,  since 
refractive  index  gradients  of  five  percent  or  better  have 
already  been  produced,  it  is  entirely  feasible  that  if  the 
precise  parabolic  change  could  be  controlled,  this  lens  could 
be  produced  today. 

Despite  the  obvious  success  of  GRIN  in  controlling  spot 
size  growth  and  image  centroid  movement,  a  penalty  in  the 
form  of  reduced  ray  intensity  has  been  paid.  That  increasing 
spot  size  performance  is  tempered  by  a  loss  of  intensity  may 
be  seen  by  comparing  "good"  GRIN  encircled  energy  plots  with 
that  of  the  HIN  lens.  Note  that  the  HIN  lens  with  N2  =  3.0 
also  loses  ~  50%  to  intensity;  see  Figure  44.  This  loss  in 
intensity  is  partially  offset,  however,  by  the  GRIN  lens 
with  an  increased  number  of  rays  transmitted  and  by  the  re¬ 
duced  3 pot  size. 

One  drawback  to  the  use  of  the  GRIN  lens  as  a  self 
sufficient,  single  element  lens  is  still  that  of  the  relative 
sizes  of  image  and  detector.  Hence,  even  though  GRIN  has 
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significantly  reduced  spot  size  below  a  =0.6,  that  size  is 

P 

still  significant  at  tilt  angles  above  0.3  radians.  In  order 
to  use  the  lens  without  a  secondary  focusing  element  or 
mirror  arrangement,  a  large  composite  sensor  array  would  be 
required. 


VII.  CONCLUSION 


Gradient  refractive  index  materials  may  be  employed  to 
design  a  pointed  seeker  lens  which  exhibits  optical  per¬ 
formance  far  superior  to  that  obtainable  with  conventional 
homogeneous  optical  material.  A  fifty  percent,  positive, 
parabolic  gradient  index  with  center  of  symmetry  interior 
to  the  lens  was  found  to  yield  the  best  performance  although 
a  five  percent  version  of  the  same  lens  was  only  very  slightly 
inferior;  this  lens  may  possibly  be  fabricated  today. 

For  objects  off-axis  by  more  them  0.3  radians  (-17.2 
degrees)  a  secondary  lens  element  may  be  required  unless  a 
large  scale  multiple  element  sensor  array  is  employed.  With 
such  arrays,  objects  off-axis  by  more  than  0.65  radians  (37.2 
degrees)  may  require  Cassegrainian  or  other  mirror  elements 
to  compensate  image  movement. 


VIII.  RECOMMENDATIONS  FOR  FUTURE  WORK 


This  thesis  investigated  a  spherically  symmetric  GRIN 
seeker  lens  with  inside  conical  surface.  Future  studies 
should  investigate: 

a)  The  elimination  of  the  singularity  problem  by 
modification  of  the  theoretical  equations  for  GRIN 
using  numerical  techniques.  See  Reference  13. 

b)  The  optimization  of  the  GRIN  lens  with  both  inside 
and  outside  curvature. 

c)  The  effect  on  lens  performance  due  to  an  attached 
shock  wave. 

d)  The  effect  on  lens  performance  when  the  object  is 
no  longer  in  the  far  field  and  wave  front  curvature 
must  be  taken  into  account. 

e)  The  performance  of  the  lens  using  wavelengths  corres¬ 
ponding  to  atmospheric  windows. 

f)  The  effect  of  a  radially  symmetric  gradient  index 
on  lens  performance. 

g)  The  feasibility  of  adding  an  anamorphic  gradient- 
index  lens  located  on  the  missile  body  to  increase 
off-axis  tracking/acquisition  capability  to  90°  and 
beyond.  See  Figure  1  of  Reference  14 . 


APPENDIX  A 


COMPUTER  LOGIC  FLOW  DIAGRAM  FOR  PROGRAM  GISL 


This  appendix  together  with  Appendices  B  and  C  describe 
the  FORTRAN  program  GISL.  GISL  may  be  used  to  design  either 
a  homogeneous  (B  =  0) ,  or  a  GRIN  lens  to  the  user's  specifi¬ 
cations  by  changing  the  design  parameter  where  indicated  by 
comments  in  the  input  section  of  the  program.  Additional 
software  required  to  run  GISL  and  plot  the  results  are  not 
described  since  these  system-dependent  procedures  do  not 
apply  elsewhere. 

The  flow  diagrams  which  follow  provide  the  reader  with 
the  information  necessary  to  follow  the  fundamental  logic 
of  the  main  program  and  subroutines  of  GISL.  Not  shown  is 
Subroutine  DIRECT  which  has  been  derived  from  Subroutine 
DIRECT  as  described  by  Amichai  [5]  with  only  slight  modifi¬ 
cations. 


.36 


* 

Here, 


the  prime  indicates  Y-aoordinate  as  transformed  into  the  grid  plane. 
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GISL  TABULAP  OUTPUT 


APPENDIX  C 


SAMPLE  TABULAR  OUTPUT  FROM  PROGRAM  GISL 
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OBJECT  PLRNE 


Shown  in  Figure  11 
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Figure  D-2.  Spot  Diagram  Corresponding  to  Object 
Plane  of  Figure  D-l 


Figure  D-3.  Object  Plan 
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Encircled  Energy  Plot  for  the  Spot 
Diagram  of  Figure  D-4 
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Figure  D-7.  Spot  Diagram  of  HIN  Lens  at  a 
Radians,  N~  =  1.5 
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Figure  D-12.  HIN  Object  Plane  at  a  -  0.4  Radians 
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Figure  D-18.  HIN  Objec 
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Figure  D-19.  Spot  Diagram  of  HIN  Lens  at  ap  *  0.6 
Radians,  N?  *  1.5 
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Figure  D-26.  Spot  Diagram  Corresponding  to  Figure  D-25 
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APPENDIX  E 


GRIN  LENS  PERFORMANCE  PLOTS  IN  THE  LOW  RANGE  OF 
- INDICES  OF  REFRACTION  (a  »  Z.Z5) 
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Figure  E-l.  GRIN  Lens  Shape  at  +50%,  OB  ■  0.05, 
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Figure  E-5.  GRIN  Lens  Shape  at  +50%,  OB  *  0.15, 
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Figure  E-6.  Grid  Plane  at  a  =0.3  for  Lens  of  Figure 
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Figure  E-7.  Spot  Diagram  for  Grid  of  Figure  E-6 
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Figure  E-9.  GRIN  Lens  Shape  at  -2%,  OB  =  0.20, 
a  s  2.25 
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Figure  E-15.  Spot  Diagram  for  Grid  of  Figure  E-14 
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E-19.  Spot  Diagram  for  Grid  of  Figure  E-18 
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Figure  E-21.  GRIN  Lens  Shape  at  +10%,  CB  =  0.20, 
a  =  2.25 
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Figure  E-23.  Spot  Diagram  for  Grid  of  Figure  E-22 
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Figure  E-25.  GRIN  Lens  Shape  at  +25%,  OB  =  0.20 
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Figure  E-27.  Spot  Diagram  for  Grid  of  Figure  E-26 
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Figure  E-31.  Spot  Diagram  for  Grid  of  Figure  E-30 
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Figure  E-33.  GRIN  Lens  Shape  at  -10%,  OB  =  0.30, 
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Figure  E-35.  Spot  Diagram  for  Grid  of  Figure  E-34 
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Figure  E-37.  GRIN  Lens  Shape  at  -25%,  OB  =  0.35, 
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Figure  E-39.  Spot  Diagram  for  Grid  of  Figure  E-38 
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Figure  E-47.  Spot  Diagram  for  Grid  of  Figure  E-46 
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Figure  E-51.  Spot  Diagram  for  Grid  of  Figure  E-50 
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Figure  E-58.  Grid  Plane  at  a  =0.3  for  Lens  of  Figure  E-57 
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Figure  E-62.  Grid  Plane  at  a  =0.3  for  Lens  of  Figure  E-61 
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Figure  E-69.  GRIN  Lens  Shape  at  -5%,  OB  =  3.50, 
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Figure  E-79.  Spot  Diagram  for  Grid  of  Figure  E-78 
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Figure  E-82.  Grid  Plane  at  a  =  0.3  for  Lens  of  Figure  E-81 
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Figure  E-83.  Spot  Diagram  for  Grid  of  Figure  E-82 
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Figure  E-87.  Spot  Diagram  for  Grid  of  Figure  E-86 
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Figure  F-43.  Spot  Diagram  for  Grid  of  Figure  F-42 
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Figure  F-51.  Spot  Diagram  for  Grid  of  Figure  F-50 
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Figure  F-54.  Grid  Plane  at  a  =0.3  for  Lens  of  Figure  F-53 
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Figure  F-67.  Spot  Diagram  for  Grid  of  Figure  F-66 
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Figure  F-80.  Encircled  Energy  of  Figure  F-79 
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Figure  F-99.  Spot  Diagram  for  Grid  of  Figure  F-98 
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Figure  F-103.  Spot  Diagram  for  Grid  of  Figure  F-102 
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Figure  F-107.  Spot  Diagram  for  Grid  of  Figure  F-106 
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Figure  F-lll.  Spot  Diagram  for  Grid  of  Figure  F-110 
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Figure  F-123.  Spot  Diagram  for  Grid  of  Figure  F-122 
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Figure  G-20.  Spot  Diagram  for  Grid  of  Figure  G-19 
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Figure  G-24 .  Encircled  Energy  of  Figure  G-23 
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Figure  H-16.  Grid  Plane  at  a  =0.5  for  Lens  of  Figure 


PLRNE 


3NU1d  133rgo-i 


FfULEO  RAYS 


PLRNE 


O.Q  0 . i  0.2  0.3  0.4  0.5 

RADIUS  FROM  IMAGE  CENTROID 


aure  h-21.  Encircled  Energy  of  Figure  H-20 


473 


OBJECT  PLflNE 


Plane  at  a  =0.7  for  Lens  of  Figure 


aNHid  i'jsrgo 


OBJECT  PLfiNE 


-IMAGE  PLANE 


H0-fll2S  167  AERODVNPHICALLV  EFFICIENT  GRADIENT  REFRACTIVE  INDEX 
HD  Hl/3  «.<■  H,SSILE  SEEKER  LENS(U)  NAVAL  POSTGRADUATE  SCHOOL 
HONTEREV  CA  H  H  CARR  OCT  82  NPS67-82-012 
UNCLASSIFIED  F/Q  16/4 


6/6 

NL 


MICROCOPY  RESOLUTION  TEST  CHART 

•  *  NATIONAL  BUREAU  OF  STANDARDS -1963 -A 


SPOT  DIAGRAM  ENERGY  DENSITY 


RADIUS  FROM  IMAGE  CENTROID 


Figure  H-27.  Encircled  Energy  of  Figure  H-26 


LIST  OF  REFERENCES 


1.  Marchand,  Erich  W.,  "Gradient- index  imaging  optics 
today",  Applied  Optics,  V.  21,  number  6,  p.  983,  15 
March  19  82 • 

2.  Kitano,  I.,  Nishigawa,  K.,  and  Momokita,  A.,  "Diffusion 
behavior  of  doped  polarizable  ions  in  glass  during  the 
ion-exchange  process".  Applied  Optics,  V.  21,  number  6, 
p.  1017-1020,  15  March  I58JT 

3.  Frazier,  Robert  L.,  Exterior  Ballastics,  Guidance  Laws 
and  Optics  for  a  Gun^Launched  Missile,  M.S.  Thesis, 

Naval  Postgraduate  School,  Monterey, California, 

December  1980. 

4.  Terrell,  James  M. ,  Conical  Lens  for  5"/54  Gun  Launched 
Missile,  M.S.  Thesis,  Naval  Postgraduate  School, 
Monterey,  California,  June  1981. 

5.  Amichai,  oded.  Sharp  Nose  Lens  Design  Using  Refractive 
Index  Gradient,  Contractor  Report.  Naval  Postgraduate 
School ,  Monterey ,  California,  June  1982. 

6.  Otten,  L.J.  and  Gilbert,  K.G.,  Aero-Optical  Phenomena, 

V.  80,  AIAA  Progress  Series  in  Astronautics  and 
Aeronautics ,  1982. 

7.  Kings lake,  Rudolf,  Lens  Design  Fundamentals,  Academic 
Press,  1978. 

8.  Hecht,  Eugene  and  Zajac,  Alfred,  Optics,  Addison-Wesley, 

1979.  - 

9.  Marchand,  Erich  W.,  Gradient  Index  Optics,  Academic 

Press,  1978.  * 

10.  Curcio,  M.E.,  "Precision-machined  optics  for  reducing 
system  complexity",  Infrared  Imaging  Systems  Technology, 
Society  of  Photo-Optical  Instrumentation  Engineers  (^PIE) 
V.  226,  p.  91-97,  April  1980. 

11.  Harney,  R.C.,  "Dual  active/passive  infrared  imaging  sys¬ 
tems",  Infrared  Imaging  Systems  Technology,  Society  of 
Photo-Optical  instrumentation  Engineers  (SPIE) ,  V.  226, 
p.  74-82,  April  1980. 


480 


12.  Robb,  Paul  N. ,  "Lena  design  using  optical  aberration 
coefficients".  Proceedings  of  1980  International  Lens 
Design  Conference  ( OS a) ,  Society  of  Photo-Optical 
Instrumentation  Engineers  (SPIE) ,  V.  237,  p.  109*118, 
1980. 


13.  Sharma,  Anurog,  Kumar,  D.A.  Vizia,  and  Ghatak,  A.K., 

"Tracing  rays  through  graded* index  media:  a  new  method". 
Allied  Optics,  V.  21,  number  6,  p.  984*987,  IS  March 


14.  Moore,  Duncan  T.,  and  Stagaman,  Joan  M. ,  "Ray  tracing  in 
anamorphic  gradient* index  media”.  Applied  Optics,  V.  21, 
number  6,  p.  999-1003,  15  March  19727 


481 


INITIAL  DISTRIBUTION  LIST 


1.  Defense  Technical  Information  Center 
Cameron  Station 

Alexandria,  Virginia  22314 

2.  Library,  Code  0142 
Naval  Postgraduate  School 
Monterey,  California  93940 

3.  Department  Chairman,  Code  67 
Department  of  Aeronautics 
Naval  Postgraduate  School 
Monterey,  California  93940 

4.  Distinguished  Professor  A.  E.  Fuhs 
Code  67Fu 

Department  of  Aeronautics 
Naval  Postgraduate  School 
Monterey,  California  93940 

5.  LTC  Rene  Larriva,  USMC 

Defense  Advanced  Research  Projects  Agency 
1400  Wilson  Blvd. 

Arlington,  Virginia  22209 

6.  Commander,  Naval  Sea  Systems  Comoand 
Naval  Sea  Systems  Command  Headquarters 
Attn:  Code  62YC 

Washington,  D.C.  20362 

7.  Mr.  Conrad  Brandts 

Naval  Surface  Weapons  Center 
Dahlgren,  Virginia  22448 

8.  Dr.  Fred  Billig 
Applied  Physics  Lab 
Johns  Hopkins  Road 
Laurel,  Maryland  20810 

9.  Commander,  U.S.  Army  Armament  Research  and 

Development  Command  (ARRAD  COM) 

Attn:  Mr.  Lou  Marino 
Dover,  N.J.  07801 

10.  Deputy  Chief  of  Staff  for  Research 
Development  and  Studies 
Headquarters  USMC 
Washington,  D.C.  20370 


No .  Copies 
2 

2 

1 

4 


2 

1 

1 

1 

1 

1 


Commander 

Development  Command 
USMC  Base 

Quant ico,  Virginia  22134 
Commander 

Artillery  Development  Command 
Fort  Sill,  Oklahoma  73503 

Captain  Herbert  M.  Carr 
2100  Raleigh  Ave. 

Austin,  TX  78703 

Dr.  Lloyd  Smith 
Code  3205 

Naval  Weapons  Center 

China  Lake,  California  93555 

Dr.  Harold  L.  Bennett 
Code  3272 

Naval  weapons  Center 

China  Lake,  California  93555 

Dr.  T.  G.  Bergman 
Code  3941 

Naval  Weapons  Center 

China  Lake,  California  93555 

Dr.  Oded  Amichai 
Ministry  of  Defense 
P.O.B.  2250 
Haifa,  ISRAEL 

Research  Administration 
Code  012a 

Naval  Postgraduate  School 
Monterey,  California  93940 

Deputy  Under  Secretary  of  the  Army 
for  Operations  Research 
Room  2E261,  Pentagon 
Washington,  D.C.  20310 

LT  James  Terrell,  USN 
Code  290 

Supervisor  of  Shipbuilding,  Conversion 
and  Repair,  USN 
Pascagoula,  MS  39567 


